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SUMMARY 

A FORTRAN computer program for the semianalytic prediction of unsteady thrust augmenting 
ejector performance has been developed, based on a theoretical analysis for ejectors. That analysis 
blends classic self-similar turbulent jet descriptions with control-volume mixing region elements. 
Division of the ejector into an inlet, diffuser, and mixing region allowed flexibility in the modelling of the 
physics for each region. In particular, the inlet and diffuser analyses are simplified by a quasi-steady- 
analysis, justified by the assumption that pressure is the forcing function in those regions. Only the 
mixing region is assumed to be dominated by viscous effects. 

The present work provides an overview of the code structure, a description of the required input 
and output data file formats, and the results for a test case. Since there are limitations to the code for 
applications outside the bounds of the test case, the user should consider TEA as a research code (not as 
a production code), designed specifically as an implementation of the proposed ejector theory. Program 
error flags are discussed, and some diagnostic routines are presented. 


CODE OVERVIEW 

TEA is a FORTRAN computer program (TEA .FOR), designed for batch execution, which reads an 
input file (TEA .DAT) to obtain the necessary program data (ejector geometry, primary nozzle conditions, 
free-stream conditions), solution control options, and file output requests for print and plot files 
(TEA. OUT and TEA.PLT, respectively). An overview of the TEA .FOR program follows; the required 
TEA .DAT input file structure is the subject of the next section. 

The main program simply serves to open the necessary I/O files and instructs the processing of 
input data, after which the governing equations are solved. 1 

A total of three subroutines are contained in the TEA.FOR file; 2 they are 

DIP - Data Input Processor 

MEC - Main Execution Control 

SPC - Single Point Calculation 


1 The solution procedure can be influenced by input file specifications. 

2 Appendix B provides a complete fisting of the TEA .FOR file. Some idea of complexity is given by 
code length: the code is 731 fines long, divided into 21 fines for the main program, 164 for the DIP, 57 
for the MEC, and 569 for the SPC. Programming style leans toward clarity rather than conciseness. 
The development platform for this code was a VAX 8800; the code has not been tested on other 
platforms. 



Figure 1 illustrates that the SPC routine is called by the MEC routine. Although the primary 
purpose of this document is to provide instructions for preparation of the TEA.DAT file (read by the 
DIP), some brief remarks on each program module are helpful in interpreting the overall program 
structure. 


Module Descriptions 

Datai inmit nrocessor . — The data input processor (DIP) is designed to obtain baseline ejector 
engineering information, free-stream properties, solution control options, and printout and plot control 
parameters from the input data file and place this information in COMMON blocks; the use of COM- 
MON blocks reduces (or eliminates) the need for READ statements in other program subroutines. 

Input data is read from TEA .DAT according to a macro command sequence. The DIP echoes 
essential program data (specified and default) so the user can conveniently identify incomplete or 
incorrect input data. In the interest of clarity, the DIP routine was not written as compactly as possible; 
this is an overriding perspective throughout the code. For ease of use, all data has a free-format 
structure — the user need only demarcate values with a comma, rather than be concerned with the 
customary Hollerith and data format (e.g., 2F12.3). 

Main mrpwitinn control . — The main execution control (MEC) routine provides the basic service of 
feeding the single point calculation (SPC) subroutine argument list with the information needed to 
compute ejector thrust at a specified point in time. Although it may appear the MEC is a trivial routine 
to include, it facilitates development of an SPC that is void of I/O for the time-dependent primary flow 
conditions. 

A previously mentioned aspect of the DIP and MEC subroutine structure is the absence of 
subroutine calls — all information is transmitted through COMMON blocks. While this may not be the 
most effective programming practice, 8 it is an approach that met the simulation specifications at the 
time the code was under development. 

Single point computation . — The single point computation (SPC) subroutine is the focal point of the 
TEA program and provides the code implementation of the transient ejector theory given by Drummond 
(ref. 1). Execution of the SPC routine is controlled by the MEC, since the SPC is structured to be a 
“drop-in” module for integration in, for instance, an aircraft engine simulation. The SPC is the largest 
program subroutine, consuming almost two-thirds of the overall TEA code (see lines 243-731 of the 
listing in Appendix B). 

Information for the SPC routine is transmitted primarily through COMMON blocks. This 
influence on the TEA computer program was the result of a requirement that the SPC argument list 
include only the primary flow conditions (m , , P 1P ) at a given time t , and return thrust and 

augmentation ratio (T , <(>) for a real-time component-level-model engine simulation. 3 4 This requirement 
is clearly simulation specific, but does provide module parameter inputs consistent with information 


3 The COMMON block chalkboard approach to programming can make program error tracking during 
diagnostics a difficult process. 

4 As discussed by Drummond (ref. 1), the motivating factor behind the development of the transient 
ejector theory was the need to develop a high-fidelity propulsion system simulation for the Advanced 
Short Take-Off and Landing (ASTOVL) development program. 
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typically available to an engine simulation code. 5 Although alternate simulation requirements may 
result in the need for a different SPC argument list, the modular code structure remains quite conve- 
nient. 

Co mm ent car ds in the SPC routine mark out blocks of calculations corresponding with calculation 
steps in the theory; an outline of the SPC solution procedure is shown in figure 2. 

Within the SPC module, there exists a decision branch between quasi-steady and transient flow 
calculations, based on the value of an execution mode switch, MODE, in the SPC argument list; this 
points to an important aspect of an ejector mixing region analysis. Since it is expected that transient 
ejector performance predictions are kicked-off from an assigned initial condition (there can be more than 
one initial condition of interest), there is a need for a built-in initialization routine for the field variables. 
A q ua si-steady solution to the governing equations is associated with a value of MODE=0; such a 
solution does not invoke time derivatives in the analysis and provides, by default, a series of steady-state 
flow computations. When MODE=l, a transient solution is invoked which reflects the implementation 
of the proposed transient ejector theory. 

System of units . — Few applications in engineering analysis appear to be exempt from the need to 
clarify the system of units involved in calculations. Typically, confusion extends from translating back 
and forth between force, mass, and pressure; the usual remark is that a consistent system of units must 
be used in analysis, and that certainly is the case here. In the present work, the Engineering English 
system of units is used. 

With the backdrop that the current work represents a research code, the use of the Engineering 
English system was appropriate for the system simulation during the code development period (in which 
this code was ultimately intended to be compatible with a propulsion system simulation that used a 
similar unit system; the Metric system was not appropriate for the simulation application). 

The Engineering English system is based on four fundamental units for length, force, time, and 
mass Units of length are in feet (ft), units of time are in seconds (s), units of force are in pounds-force 
(lbp, and units of mass are in pounds-mass (lb m ). Compatibility of units is provided by the introduction 
of the gravitational constant g c = 32.174 into Newton's law (gj, is numerically, but not dimensionally 
equal to the gravitational constant). Although the work of Drummond (ref. 1) provides extensive 
examples on the appropriate units for each variable involved in the analysis, the following summary is 
relevant to input file preparation: 


Description 

Mnemonic 

Value = Units 

Temperature 

T 

= °R 

Pressure 

P 

= lb/ft 2 

Velocity 

V 

= ft/s 

Density 

RHO 

= lb m /ft 3 

Mass flowrate 

MDOT 

= lb m /s 

Specific heat ratio 

GAMMA 

1.4 = non-dim 

Gas constant 

RBAR 

1715.66 = ft 2 /s 2 R 

Angle of attack 

ALPHA 

= degrees 


5 This leads to ce rtain calculations that normal ejector codes would not have a requirement for, or the 
use of (e.g., a total pressure versus static pressure specification at the secondary inlet). 
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In principle, there is nothing to prevent the user from using, say, the Absolute MKS system of units 

(g has a value of unity) since equations involving the value of gc do so symbolically. Illustrations of this 

are shown in lines 118 and 122 of the SPC routine - g c is set at the top of the SPC module (line 67) and 
could be reassigned as desired. 


INPUT DATA FILE PREPARATION 
Typical File Structure 

Data processing for the TEA. FOR program involves extensive use of macro commands; macros are 
keywords placed in the input file which permit a flexible input data file structure. Appendix A provides 
an example 6 input file structure. 

Although the macro command structure can be perceived as a somewhat clumsy programming 
practice, the result is a very lucid input file structure. Another benefit to the macro approach is that the 
blocks of data are not required to be placed in any particular order, and in some cases they need not 
appear at all (as long as the program defaults are acceptable to the user). A subtle feature of the macro 
command approach is that the program will ignore any lines that do not begin with a recognized 
keyword — as such, the user can include as many descriptive comment cards as desired at the top of the 
file for a full description of the case study at hand. 7 


Macro Commands and Options 


Wi thin the data input processor (DIP), solution instructions, data input, and data output 
specifications are read from the TEA. DAT file through a macro command structure. As mentioned 
previously COMMON statements are extensively used as placeholders in memory to shuttle data from 
one routine to another. Except for the SPC, none of the subroutines have argument lists since the 
necessary information is provided by the COMMON statements. For this type of data transfer, sub- 
program macros are very useful. 

Data to be processed within the DIP block is subdivided according to a logical characterization of 
the data to be handled. Currently, eight keywords are used in the subdivision of data: 


■BACK' PRESSURE 
'CASE* 

'COEF ' FICIENT 
•DISP’LAY 
TREE 'STREAM 
'GEOM'ETRY 
TRIM' ARY 
'STOP' 


Back pressure specification for ejector discharge (optional) 
Case banner for printed output (optional) 

Empirical coefficients for mixing and entrainment 
Display control for printed and plotted output (optional) 


Free-stream fluid properties 

Geometry of the ejector 

Primary nozzle discharge conditions 

Stop command for data input processing (optional) 


6 This file corresponds with the test case to be discussed later, for now, however, it is the file structure 
rather than the numbers themselves that are of interest 

7 After several weeks pass, it is easy to forget why a particular case was of interest. 
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Quotes around the first four letters of each keyword are used to indicate the minimum portion of the 
word to be included for data file preparation. Only four of the eight data blocks are essential to include 
in the input file; optional sets are marked accordingly. Except for STOP, which (if used) must be the last 
keyword encountered in the DIP data block, the keywords and associated data can occur in any se- 
quence desired. Descriptions of the input parameters associated with each macro command are given 
below (in alphabetical order). 

BACK: Ejector back-pressure (optional) . — Aircraft installation aerodynamics often dictate that 
the back pressure of an installed ejector may be different than atmospheric; to accommodate this, the 
BACK block allows the prescription of a nonatmospheric exit pressure. 


Line 

Variable 

Description 

1 

2 

’BACK 

’PBACK 

Keyword for back pressure input block 

Value for ejector back pressure [Default: PBACK=PINF] 


This data block is optional if the back pressure is equal to that of the ambient since the program default 
is to assign the backpressure to the value of the free-stream static pressure (as specified and discussed 
later in the FREE stream block). PBACK has units of lb f /ft. 2 

CASE: Banner for printed output (optional) . — The descriptive banner is the header for printed 
output This data block is composed of two lines: 


Line 

Variable 

Description 

i 

'CASE' 

Keyword for banner input block 

2 

’RUNID' 

Banner for printed output [Default: RUNID='TEA'] 


The CASE block is optional; the program provides a default banner of TEA if this data block is omitted. 

COEF: Coefficients for solution ty pe- mixing, kinetic energy entrainment — As discussed in the 
work of Drummond (refs.l and 2) the ejector theory implemented in this work employs two empirical 
constants, the flow mixin g effectiveness, BETA (B), and the coefficient for the kinetic energy transfer 
function, CSIGMA (C s ). 


Line 

Variable 

Description 

i 

’COEF' 

Keyword for coefficient data block 

2 

'TITLE' 

Comment fine with descriptive titles for variables on next line 

3 

'ITYPE', 

Type of solution desired (0 = Steady; 1 = Transient) 


'BETA', 

Mixing effectiveness factor (1.0 < BETA < 2.0) 


'CSIGMA' 

Kinetic energy coefficient (0.0 < CSIGMA < 1.0) 


In the interest of making direct comparisons between quasi-steady and full transient flows, this block 
introduces the option of enforcing a quasi-steady flow assumption. What essentially happens is that the 
steady-flow initialization routine is invoked for all point calculations emanating from the MEC. The 
default integer has a value of ITYPE=1 and marks the full transient flow as the execution mode; a value 
of 0 eliminates time derivatives to provide the quasi-steady option (ITYPE sets the value of the MODE 
parameter discussed previously in the SPC overview). 
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TYTSP: Printout of intermediate computations (optional) . — For detailed investigations of transient 
performance calculation results, three parameters can be set to control printed output and plot file 
information. 


Line 

Variable 

Description 

1 

'DISF 

Keyword for display data block 

2 

'TITLE' 

Comment line with descriptive titles for variables on next line 

3 

'ICHECK, 

Integer for printed calculation output [Default: ICHECK=100] 


'INDEX', 

Indexing for printout loops [Default: INDEX=1] 


’IPLOT 

Index for plotfile data [Default: IPLOT=0] 


Values of ICHECK can range between 0 and 1000, where 0 provides no printed output at all, 50 provides 
only a print summary from the DIP and MEC routines, and a value of 100 (the default) provides detailed 
lis ting of SPC results as each time step calculation is completed. 8 For diagnostic applications, 

ICHECK = 500 provides detailed velocity initialization data, and at ICHECK = 1000, self-similar profiles 
and secondary velocities are printed for every iteration within each time step. The INDEX parameter is 
useful if the program has operated on, for instance, 500 primary nozzle states, and it is desired to only 
printout every 10th computed result; INDEX must have a value greater than or equal to 1. The IPLOT 
parameter dictates whether a plot file is to be created, and the intervals at which data is written to the 
plot file. For the default value of IPLOT = 0, no data is stored in the TEA .PLT file. For IPLOT = 1, 
primary flowrate, thrust, and thrust augmentation are stored at each time interval; for values of 
IPLOT = N, every Nth data point is stored (much like the INDEX function). 

Plot file output parameters are controlled by the WRITE statement at lines 231-235 in the code. At 
present, only parameters passed through the SPC argument list are amenable to place in the plot file; 
other parameters could be used if the appropriate WRITE statements were written to replace the 
listing segment of code. 

FREE: FVe^gt-rpam fluid properties (optional) . — Characteristic (free-stream) flow properties and 
thermodynamic data are contained in the FREE stream data block. 


line 

Variable 

Description 

i 

’FREE’ 

Keyword for free-stream data block 


2 

'TITLE' 

Comment line with descriptive titles for variables on next line j 

3 

TINT", 

Free-stream static temperature 

[= 530 n; default value] 


PINT, 

Free-stream static pressure 

[= 2116.8 lb/ft 2 ] 


'GAMMA', 

Specific heat ratio 

[= 1.4 nondim] 


'RBAR', 

Universal gas constant 

[= 1715.66 ft 2 /s 2 R] 


UINF, 

Free-stream velocity 

[=0.0 ft/s] 


'ALPHA' 

Angle of attack 

[= 90.0 deg] 


To simplify the implementation of the ejector theory, the primary and secondary ejector streams are 
assumed to have the same gas composition (this is not a restriction of the theory). 


8 The rationale for the ICHECK values of 0, 50, 100 (instead of 1, 2, 3) extends from the prototype code 
effort to interface with the airframe simulation. 
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GEOM: Elector geometric characteristics . — Geometric approximations for the ejector are con- 
tained in the GEOM data block. All numeric data must be proceeded by a title. Data is expected in the 
following order 


Line 

Variable 

Description 

1 

'GEOM' 

Keyword geometric data input block 

2 

'TITLE' 

Comment line with descriptive titles for variables on next line 

3 

■X0,Y0' 

Inlet (station 0) x- and y-dimensions 

4 

TITLE’ 

Comment line with descriptive titles for variables on next line 

5 

XI, Yl, 

Primary nozzle (Stn 1) x- and y-dimensions x-section 


NX1,NYT 

Number of nozzles in the x-,y- directions 

6 

TITLE' 

Comment line with descriptive titles for variables on next line 

7 

’X2.Y2’ 

Secondary flow (station 2) x- and y-dimensions 

8 

'TITLE' 

Comment line with descriptive titles for variables on next line 

9 

ra/Y^ 

Mixing region (station 3) x- and y-dimensions 

10 

TITLE' 

Comment line with descriptive titles for variables on next line 

11 

■X4,Y4’ 

Diffuser exit (station 4) x- and y-dimensions 

12 

TITLE' 

Comment line with descriptive titles for variables on next line 

13 

'zi,Nzr 

Total mixing region length and number of finite volumes 


Figure 3 illustrates the expected geometric layout of the ejector, specific data was shown previously in 
Appendix A. 

PRIM- Primary nnzsle discharge conditions . — This block specifies the primary nozzle discharge 
velocities and thermodynamic properties. 


line 

Variable 

Description 

i 

'PRIM 

Keyword for primary nozzle data block 

2 

TYPE' 

Keyword for nature of data to follow: 



‘SPEC 7 -> Data specified for eveiy t 



'STEP -> Step function specification 


If TYPE = 'SPEC', then the following input data sequence is expected: 


Line 

Variable 

Description 

3 

'TITLE' 

Comment line with descriptive titles for variables on next line 

4,...N 

’T(I), 

Time stamp 


DMPDTO), 

Primary flowrate 


T1P0(D, 

Stagnation temperature 


P1P(I)' 

Static pressure 

N+l 

0 

Numbers) or characters) to trigger end of record 


Note that there are four data values on each of the 4,...N lines to simulate the primary nozzle conditions 
as a function of time. The N+l reading must be a number or letter format that will terminate input by 
triggering an error in the READ loop (a simple single value of‘0’ works fine; see Appendix A). 
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If TYPE = 'STEF, then the following input data sequence is expected: 


Line 

Variable 

Description 

3 

TTTT.F.' 

Descriptive comment for variable list 

4 

TX1) 

Starting time 


DMPDTU) 

Starting primary flowrate 


T1PCK1) 

Starting stagnation temperature 


P1P(1)' 

Starting static pressure 

5 

TITLE' 

Comment line with descriptive titles for variables on next fine 

6 

•Z1TIME' 

Time at which step-change occurs 


’Z1DMDT 

New primary flowrate 


'ZITIPO' 

New stagnation temperature 


'ZIP IF 

New static pressure 

7 

TITLE' 

Comment fine with descriptive titles for variables on next line 

8 

TSTEP.N 1 

Time-step size and number of time steps to take 


QTnp Statement (nntional) . — The end of input data is marked by the presence of the QUIT card. 


Line 

Variable 

Description 

i 

'STOP 

Keyword to exit DIP routine 


ILLUSTRATIVE CASE 
Background 

Research on design methodologies for integrated aircraft and flight control systems requires 
accurate subsystem component simulations. In principle, these simulations must mimic steady-state 
and transient component effects. Of particular interest in recent years is the development of real-time 
simulations for short take-off vertical landing (STOVL) aircraft.* The thrust augmenting ejector shown 
in figure 4 is considered a potentially valuable propulsion subsystem element for the powered-lift aspect 
of STOVL aircraft. 

Steady-flow ejector performance predictions are typically the centerpiece of an ejector design cycle, 
and transient flow behavior is most often described through a quasi-steady flow approximation. The 
latter essentially assumes away the transient behavior issue since the conjecture is the ejector will 
instantaneously provide whatever thrust is demanded. 

To establish the validity of the quasi-steady flow assumption, one must be able to quantify the 
elapsed time for the ejector to adapt to a new thrust level, and establish that the ejector delay is less 
than the control system update period. In air-breathing propulsion systems, this assessment can take 
the form of a comparison between the ejector frequency response and that of the propulsion control 
system; see figure 5. The purpose of the TEA program is to assist in the assessment of the ejector 
response rate; the specific case study used to examine the STOVL response rate is presented next. 


» For more extensive details on the simulation of interest, see Drummond (ref. 2), or Drummond and 
Ouzts (ref. 5). 
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Characteristic Test Case 


In the absence of appropriate 10 data from transient flow ejector experiments or from advanced time- 
dependent CFD calculations, verification that the proposed ejector analysis can provide reasonable 
thrust predictions must be left to engineering judgement. Because of this, a familiar forcing function 
must be used. In the present work, the system response to a step-function input is not only a 
characteristic transient case study, but the scenario is coincident with typical STOVL ejector 
application. For demonstration purposes the ejector system response to a step-change in primary nozzle 
flowrate from 18.7 to 21.85 lb n /sec is chosen because (1) experimental steady-state data at each of these 
operating points is available, mid (2) the 17 percent change in primary flowrate is beyond the small- 
perturbation range; this exercises the system nonlinearities. Two versions of the TEA.DAT input file 
were prepared 11 for this case study and are given in Appendix C. 

Geometric approximations for the ejector inlet, mixing region, and diffuser are presented in 
table I note the 3-D to 2-D conversion required for the primary nozzle area. This conversion in nozzle 
geometry is required so that the input data will be consistent with the mathematical formulation of the 
ejector mixing layer. 


Results 

For the mixing region, a finite-volume length of 0.18 ft and a characteristic mixing region velocity 
of 500 ft/s, the characteristic time step for computations is approximately 0.4 ms. To avoid infringing on 
this stability limi t a computational time step of 0.1 ms was chosen; 100 time steps provided the 
necessary interval for examination of the step-function test case. 

A de tail ed listing of the program printout is given in Appendix D, though the results for each of the 
100 time steps are not shown (for brevity). 

The empirical coefficient in the transient analysis (required for calibration of the primaiy-to- 
secondary kinetic energy exchange mechanism), was selected to match the asymptotic thrust prediction 
with the quasi-steady value at the new set point. Figure 6 illustrates the predicted ejector thrust profile 
for a coefficient C o of approximately 0.25. It appears the 2 millisecond residence time of the flow 
(elapsed time for the primary nozzle flow to reach the diffuser exit plane) is slightly less than the 3 
millis econd interval for the thrust to reach a new maximum. Oscillations in thrust after that point 
appear to settle out in about 5 milliseconds. 

An (initially) unexpected feature of the thrust profile is the dip in thrust immediately following the 
step-change in pr imar y nozzle efflux. Examination of the field variable profiles reveals this is not a 
numerical problem, but that the increase in static pressure associated with the instantaneous change m 
driving flow temporarily impedes the secondary flow. After a short period, the secondary flow kinetic 
energy builds to overcome this effect and then continues in the intuitively expected manner. 


10 A test was conducted by Drummond (ref. 3). specifically to collect transient ejector data, but the 
forcing function had a characteristic time much greater than the time scale of interest in simulation 
validation. 

11 The file identifier TEA .DAT is hard-coded into the program for simplicity in program development. 
As such, the examples given in Appendix C are titled CASE1.DAT and CASE2.DAT, for which it is 
assumed the user will rename the case of interest to be TEA .DAT, prior to execution. The user can 
automate this process through a simple command language script. 
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Discussion 


Steady-state ejector performance is well-known to be highly sensitive to design conditions, so it 
should come as no surprise that the transient performance prediction task is as much an art as it is a 
science. Some aspects of the illustrative test case runs are useful to discuss. 

Entrainment coefficient variation. — A critical step in the implementation of the transient ejector 
theory is the identification of the appropriate value of the empirical entrainment coefficient C,. For a 
known ejector step-function response, the desired C. value is the one that provides a match between the 
predicted asymptote and the known thrust at the second state. As an example, thrust predictions for 
various C values are shown in figure 8; this chart illustrates the idea that the TEA code was not 
designed to provide the optimal prediction of steady-state ejector performance 12 — rather, that the code 
focuses on analysis of transient flow characteristics. (Figure 7 differs slightly from figure 6 m thrust 
prediction due to a slight modification in the implementation of the theory, figure 7 is the more recent 

result). 

Rack pressure effects . — It is possible for several physical reasons that the ejector backpressure 
will differ from the (upstream) ambient as, for instance, in the case of installation effects or curved 
mixing regions (asymmetric mixing region). To this end, a performance feature to explore is 
backpressure limits. The principal goal of increasing the backpressure for any parametric investigation 
is to establish the threshold at which the mixing region discharge energy is inadequate to overcome the 
back pressure — in that situation the program notifies the user that a fatal error has occurred (see 
Appendix A) For the illustrative test case, the backpressure cannot exceed 5 percent over ambient 
Error codes are discussed further in the ERROR CODES AND EXECUTION PROBLEMS section. 

Performance sensitivity to nozzle geometry . — Ejector performance is very sensitive to ejector 
geometry. In particular, performance is most readily influenced by the primary nozzle discharge area. 
For the baseline case, a 20-percent increase in the primary nozzle area yields a 14-percent decrease m 
primary nozzle discharge velocity, but a 20-percent decrease in the primary nozzle area causes some 
difficulty in initialization of the program. As illustrated in Appendix A, the initialization difficulty is 
presented to the user through an error message for a negative entrained velocity. However, during the 
initialization process, a simple fix is to reset the entrained kinetic energy to a small positive number. 
This patch seems to work well in the few cases it has been invoked. 

A point to be made here is that the ejector transient flow solution involves a highly complex 
solution logic, and that to encounter an error is not so much a reflection that the input data is incorrect 
but that an idiosyncrasy with the solution method may have been encountered. 

The discussion in the next section suggests that previous problems encountered with execution 
have been fixed by altering input data values. Recommendations for algorithm changes or 
improvements to the theory are outside the scope of this work. 


12 B u t a sample routine is included to initialize the transient flow calculation procedure. 
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ERROR CODES AND EXECUTION PROBLEMS 


Error Codes 

At several points in the previous discussion, it has been noted that program execution may 
terminate under certain conditions. Five types of errors are currently flagged by the computer, although 
one could easily inspect the code and potentially identify many more diagnostic flags, those included 
presently extend from execution experiences with the test case described earlier. Error codes, the current 
flag location (given by the SPC line number), and a brief description of the problem are given below: 


Error 

code 

SPC 

line number 

Description of problem 

001 

272 

Imaginary value (B 2 - 4AC < 0 in the quadratic equation) has been 
computed for the centerline velocity during initialization of the field 
variables. If ICHECK has been set to 500, then the XCHECK2 routine 
will assist in the diagnostics by printing a listing of the variables 
related to the calculation. 

002 

416 

A negative value for the nominal velocity V* has been computed (V* is 
equal to the entrained velocity V e cubed). In this case, rates of change 
in KE have diminished entrainment rates to unrealistically low 
(negative) levels. The XCHECK3 routine is automatically invoked as a 
first pass at fixing the problem. (This was discussed in the test case 
description; see Performance Sensitivity to Nozzle Geometry). 

003 

448 

Unrealistic value for density, the value is less than zero. 

004 

449 

Unrealistic value for density, the value is greater than one. 

005 

452 

Unrealistic value for temperature, the value is less than zero. 

006 

466 

Stagnation pressure at station 4 is less than the static pressure; this 
results in an unrealistic discharge exit velocity. Either the ejector back 
pressure is too high or the inlet Mach number is too low. 


These problems are generally symptomatic of the divergence of the transient solution. 


Execution Problems 

Unrealistic or divergent transient flow solutions may be the result of one of the following: 

(1) An incorrect starting solution (check input data file), 

(2) Too large a time step (modify the input data time history), or 

(3) An inappropriate expression for the function for exchange of kinetic energy between the 
p rimar y and secondary fluid streams. 
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The last item is an important one since the entrainment function relies heavily on the assigned value of 
CSIGMA (see line 185 in SPC). More details on this function and its origin are described in Drummond 

(refs. 2 and 4). 


Optional Debug Routines 

Appendix B includes four routines that may be useful for printout of intermediate information 
during the debugging process. These routines are invoked through INCLUDE statements (to simplify 
their omission or modification, if desired); typically, these diagnostic routines are unnecessary if the 
program executes without error. 


Routine 

name 

SPC 

line number 

Description 

XCHECKl 

203 

Printout of self-similar profile coefficient values 

XCHECK2 

273 

Printout of data involved in the centerline velocity calculation; 
used in conjunction with SPC error code number 1. 

XCHECK3 

417 

Corrective procedure for negative V* result; used in conjunction 
with SPC error code number 2. 

XCHECK4 

494 

Entrained velocity error and discharge mass error printouts, 
used to monitor convergence of a solution at a given time step. 


CONCLUSION 

Understanding and quantifying ejector fluid-dynamic characteristics is central to design decisions 
involved in propulsion system applications. The present work provides instructions for a FORTRAN 
computer program implementation of a semianalytic prediction of unsteady thrust augmenting ejector 
performance developed by Drummond (ref. 1). That analysis blends classic self-similar turbulent jet 
descriptions with control-volume mixing region elements. Division of the ejector into an inlet, diffuser, 
and mining region allowed flexibility in the modelling of the physics for each region. 

An overview of the code structure, a description of the required input and output data file formats, 
and the results for a test case were presented. It is emphasized that there are limitations to the code for 
applications outside the bounds of the test case, so the user should consider TEA as a research code 
designed specifically as a test of the proposed ejector theory. For program applications outside the test 
case bounds, some error flags and diagnostic routines are presented and discussed. 
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APPENDIX A - CODE SAMPLES 


TEA.DAT: SAMPLE 

C CASE1.DAT — 

C 

C Case study input file for TEAFOR transient ejector analysis 
C routine. Here a single step function is prescribed that 
C corresponds with one of the DeHaviland/PLF tests run in June 1987. 

C This file assumes that the default BACK and DISP parameters in the TEAFOR 

C source code are appropriate. 

C 


CASE banner. 

Test #2, Step change in primary flow for 87 DHC ejector design. 


GEOMetric approximation for the ejector. 


XO 

0009.36000 

XI 

0000.10801, 

X2 

0000.78000, 

X3 

0009.36000, 

X4 

0009.36000, 

Z1 

0000.90000, 


YO 

0001.25000 

Y1 

0000.10801, 

Y2 

0000.34700 

Y3 

0001.04100 

Y4 

0001.87500 

NZ 

00005 


NX1 NY1 
00012, 00003 


COEFficients for solution execution. 
1TYPE BETA CSIGMA 

00001, 0001.0550, 0000.30000 


FREE stream fluid conditions. 

TINF PINF GAMMA RBAR UN IF ALPHA 

0542.80000, 2067.84000, 0001.40000, 1714.54130, 0000.00000, 0090.00000 


PRIMary nozzle discharge conditions. 
SPECify states at each time step. 


T 

0000.00010, 

0000.00020, 

0000.00030, 

0000.00040, 

0000.00050, 

0000.00060, 

0000.00070, 

0000.00080, 

0 


MDOT1P 

0018.70000, 

0018.70000, 

0018.70000, 

0018.70000, 

0018.70000, 

0021.85000, 

0021.85000, 

0021.85000, 


T1P0 

0769.70000, 

0769.70000, 

0769.70000, 

0769.70000, 

0769.70000, 

0760.60000, 

0760.60000, 

0760.60000, 


PIP 

2219.78000 

2219.78000 

2219.78000 

2219.78000 

2219.78000 

2250.92000 

2250.92000 

2250.92000 


C 
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ARIEL$ RUN TEA 

$$$ TEA FATAL ERROR# 6, $$$ 

%MTH-F-SQUROONEG, square root of negative value 
user PC 000242C1 

%TRACE-F-TRACEBACK, symbolic stack dump follows 


module name 

routine name 

line 

SPC 

SPC 

396 

MEC 

MEC 

38 

TEA 

TEA 

16 


ARIEL$ 


ARIEL$ RUN TEA 

$$$ TEA FATAL ERROR # 2, $$$ 


*** ERROR DURING CALCULATIONS FOR T = 

DKE = 976309.563 

DKEOLD = 2096179.500 
SIGMA = 71.929 

PLOSS = 11389.766 

ENERGY1S = 1076481.875 
CHANGE = -1119870.000 
VSTAR = -1261377.625 

*** PROCEED W/ CHANGE = -0.999ENERGY1S 

*** NEW VSTAR IS = 31295.961 


FORTRAN STOP 
ARIEL$ 


0.00040 
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APPENDIX B - TEA PROGAM LISTING 

00001 C TEA.FOR: Transient Ejector Analysis -Ver.2.1 15-AUG-91 

00002 0============================================================== 

00003 C Revised: 15-Feb-93 

00004 C Colin K. Drummond 

00005 C Propulsion Systems Division 

00006 C NASA Lewis Research Center 

00007 C Ref: “A Control-Volume Method for Analysis of Unsteady Thrust 

00008 C Augmenting Ejector Flows,” NASA CR- 182203. 

00009 0============================================================== 

00010 PROGRAM TEA 

00011 C 

00012 OPEN (UNIT= 5, FILE= 'TE AD AT', STATU S= OLD’) ! Input data file 

00013 OPEN (UNIT= 6, FILE= ’TEA.OUT’,STATUS= ’NEW') ! Output data file 

00014 OPEN (UNIT=10, FILE= 'TEAPLT', STATUS= 'NEW') ! Plot data file 

00015 CALL DIP ! Data Input Processor 

00016 CALL MEC ! Main Execution Control 

00017 CLOSE (UNIT= 5) 

00018 CLOSE (UNIT= 6) 

00019 CLOSE (UNIT=10) 

00020 STOP 

00021 END 
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00001 Csss*t«K*Mt«a* = = = : = s=:=="" 5S = ss = ss = = sss8SSS8S5S8S,sss,tsssssssS2 

00002 SUBROUTINE DIP 

00003 C 

00004 C Data Input Processor (DIP) 

00005 C 

00006 PARAMETER (M=5000) 

00007 CHARACTER*80 RUNID 

00008 CHARACTERM TITLE(20) 

00009 DIMENSION T (M) , DMPDT (M) ,T1P0(M) ,P1P(M) 

00010 COMMON /CHECK/ I CHECK, INDEX 

00011 COMMON /COEFX/ BETA.CSIGMA 

00012 COMMON /ITYPE/ ITYPE 

00013 COMMON /PLOTS/ I PLOT 

00014 COMMON /LIST1/ XO, YO, 

00015 4 XI, Yl, NX1, NY1, 

00016 4 X2, Y2, 

00017 4 X3, Y3, 

00018 4 X4, Y4, 

00019 4 Zl, NZ1, DZ 

00020 COMMON /LIST2/ GAMMA, RBAR, TINF, PINF, PBACK 

00021 COMMON /LIST3/ ALPHA, UINF 

00022 COMMON /FL0W1/ ISW, N, T, DMPDT, T1P0, PIP 

00023 COMMON /FL0W2/ Z1TIME, Z1DMDT, Z1T1P0, Z1P1P, TSTEP 

00024 C 

00025 C DEFAULT DATA 

00026 ALPHA = 90.00 1 Angle of attack of aircraft 

00027 GAMMA = 1.40 ! Specific heat ratio 

00028 ICHECK = 100 ! Print detailed time step output 

00029 INDEX = 1 ! Print every 1*1 times 

00030 I PLOT =0 ! No plotted output 

00031 ITYPE = 1 ! Full transient solution (= 0, steady-state) 

00032 PINF * 2067.84 ! Free-stream pressure, lbf/ft2 (14.7 psla) 

00033 PBACK = PINF 

00034 RBAR = 1714.54 1 Gas constant for AIR, ft2/s2R 

00035 RUNID = 'TEA' ! Default RUNID 

00036 TINF « 542.80 ! Free-stream temperature, R (70degF) 

00037 UINF = 0.0 ! Free-stream velocity, ft/s (stationary ejector) 

00038 C 

00039 C READ INPUT DECK 

00040 10 READ (5, ' (20A4) ' ,END=200) TITLE 

00041 IF (TITLE(l).EQ.'BACK') GO TO 100 

00042 IF (TITLE(l).EQ. 'CASE') GO TO 110 

00043 IF (TITLE(l) .EQ. ' COEF ' ) GO TO 120 

00044 IF (TITLE(l).EQ.'DISP') GO TO 130 

00045 IF (TITLE(l) .EQ. ' FREE* ) GO TO 140 

00046 IF (TITLE(l).EQ.'GEOM') GO TO 150 

00047 IF (TITLE(l) .EQ. 'PRIM') GO TO 160 

00048 IF (TITLE(l).EQ. 'STOP') GO TO 200 

00049 GO TO 10 

00050 C 

00051 C BACK pressure path to read PBACK specification. 

00052 C *** Note: This Is ONLY needed If PBACK.NE.PINF 

00053 100 READ (5,*) PBACK 

00054 GO TO 10 

00055 C 

00056 C CASE path to read RUNID banner. 

00057 110 READ (5,' (A)') RUNID 
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00058 GO TO 10 

00059 C 

00060 C COEFflclents path to read solution limits and coefficients 

00061 120 READ (5,*) 

00062 READ (5,*) ITYPE,BETA,CSIGMA 

00063 GO TO 10 

00064 C 

00065 C DISPlay path for printout and plotting options. 

00066 130 READ (5,*) 

00067 READ (5,*) ICHECK, INDEX, I PLOT 

00068 GO TO 10 

00069 C 

00070 C FREE stream path to read flow properties. 

00071 140 READ (5,*) 

00072 READ (5,*) TINF.PINF, GAMMA, RBAR.UINF, ALPHA 

00073 PBACK = PINF 

00074 GO TO 10 

00075 C 

00076 C GEOMetry path to read baseline ejector geometry. 

00077 150 READ (5,*) 

00078 READ (5,*) XO.YO 

00079 READ (5,*) 

00080 READ (5,*) X1,Y1,NX1,NY1 

00081 READ (5,*) 

00082 READ (5,*) X2,Y2 

00083 READ (5,*) 

00084 READ (5,*) X3.Y3 

00085 READ (5,*) 

00086 READ (5,*) X4,Y4 

00087 READ (5,*) 

00088 READ (5,*) Z1,NZ1 

00089 GO TO 10 

00090 C 

00091 C PRIMary path for specification of primary nozzle state. 

00092 160 READ (5, * (20A4) ' ) TITLE 

00093 IF (TITLE(l) . EQ. 'SPEC') GO TO 161 

00094 IF (TITLE(l) . EQ. 'STEP') GO TO 164 

00095 C — Read step-by-step nozzle condition matrix — 

00096 161 ISW = 0 

00097 READ (5,*) 

00098 DO 162 1=1 ,M 

00099 READ (5,*,ERR=163) T(I) ,DMPDT(I) ,T1P0(I) ,P1P(I) 

00100 162 CONTINUE 

00101 163 N = 1-1 

00102 K = N 

00103 GO TO 10 

00104 C — Read parameters for generating a step function — 

00105 164 ISW = 1 

00106 READ (5,*) 

00107 READ (5,*) T(l) ,DMPDT (1) , T1P0(1) , P1P(1) 

00108 READ (5,*) Z1TIME, Z1DMDT, Z1T1P0, Z1P1P 

00109 READ (5,*) 

00110 READ (5,*) TSTEP.N 

00111 K = 1 

00112 GO TO 10 

00113 C 

00114 C STOP reading data and print data summary (If ICHECK. GE. 50) 
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00115 

00116 

00117 

00118 

00119 

00120 
00121 
00122 

00123 

00124 

00125 

00126 

00127 

00128 

00129 

00130 

00131 

00132 

00133 

00134 

00135 

00136 

00137 

00138 

00139 

00140 

00141 

00142 

00143 

00144 

00145 

00146 

00147 

00148 

00149 

00150 

00151 

00152 

00153 

00154 

00155 

00156 

00157 

00158 

00159 

00160 
00161 
00162 

00163 

00164 


200 IF (ICHECK.LT. 50) GO TO 999 

WRITE (6,700) RUNID, ICHECK,X0,Y0,X1 , Y1,NX1,NY1,X2,Y2, 

& X3,Y3,X4, Y4,Z1,NZ1, BETA, CSIGHA, TINF, PINF, 

& GAMMA, RBAR.UINF, ALPHA, PBACK 

IF (ISW.EQ.O) THEN 
WRITE (6,706) 

WRITE (6, '((7X,I3,4(2X,F10.5))) ’) 

& ((I,T(I),DMPDT(I),T1P0(I),P1P(I)),I=1,K) 

ENDIF 

IF (ISW.EQ.l) 

& WRITE (6,705) T(1),DMPDT(1),T1P0(1),P1P(1), 

& Z1TIME,Z1DMDT,Z1T1P0,Z1P1P,TSTEP 

WRITE (6,703) N 
IF (ITYPE.EQ.O) WRITE (6,704) 

GO TO 999 


700 FORMAT (//,20X, '***** TEA INPUT *****',///3X,A, 

& //,5X,' CHECK PRINT LEVEL IS ',13, 

& //, 5X, 'EJECTOR GEOMETRIC CHARACTERISTICS, ’ 

& /,7X, 'XO * 1 ,F10.5,5X, 'YO = \F10.5, 

& /,7X,'X1 = ' ,F10.5,5X, *Y1 = ',F10.5, 

& / , 7X , 'NX1= ',I5,10X, 'NY1= ',15, 

& /,7X,'X2 = ' ,FI0.5,5X, ' Y2 = \F10.5, 

& /,7X,'X3 = ' , F10.5.5X, 'Y3 = ',F10.5, 

& / , 7X , 'X4 = * ,FI0.5,5X, 1 Y4 ■ ',F10.5, 

& /,7X,'Z1 = ' ,F10.5,5X, 'NZ = ',I5,5X, 

& //,5X, 'MIXING LOSS CORRECTION FACTOR' , 

& /,7X, 'BETA = ' , F10.5, 

& //, 5X, 'KINETIC ENERGY MIXING COEFFICIENT', 

& /,7X, 'CSIGMA * ' ,F10.5 

& // , 5X, 'FREE-STREAM THERMODYNAMIC PROPERTIES', 

& /,7X, 'STATIC TEMPERATURE = ',F10.5,' DEG-RANKINE' , 

& /,7X, 'STATIC PRESSURE = ',F10.5,' LB-F/FT3', 

& /,7X, 'SPECIFIC HEAT RATIO = '.F10.5, 

& /,7X, 'GAS CONSTANT, AIR = '.F10.5, ' FT2/S2-R', 

& /,7X, 'FREESTREAM VELOCITY = '.F10.5,' FT/S', 

& /,7X, 'ANGLE-OF-ATTACK = '.F10.5,' DEGREES', 

& //,5X,'EJE BACK PRESSURE IS ',F10.5,' LBF/FT2' , 

& //,5X, 'PRIMARY NOZZLE DATA',/) 

703 FORMAT (//5X, 'PROGRAM WILL OPERATE ON ',15,' DATA SETS') 

704 F0RMAT(//5X, 'QUASI-STEADY SOLUTION INVOKED') 

705 FORMAT ( 7X,'STEP FUNCTION FOR PRIMARY NOZZLE DISCHARGE', 

& /,7X, 1 T 1 ,4X, ' DMPDT* ,7X, 'T1P0 ’,5X,' P1P',4X, 

& /,6X, ' (sec) ' ,4X, ' (Ibm/s) ' ,5X, ' (degR) • ,4X, ' (lbf/ft2) ’ , 

& /,4(2X,F10.5) ,/,4(2X,F10.5),//,9X, 'TSTEP ■ \F10.5) 

706 FORMAT ( 

&7X, 'DATA' ,7X, ' T ',4X,' DMPDT* ,7X, 'T1P0 ',5X,' P1P',4X,/, 

&7X, 'SET# ' ,6X, ' (sec) ' ,4X, ' (lbra/s) * ,5X, • (degR) * ,4X, • (lbf/ft2) ') 
999 RETURN 
END 
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00001 

00002 

00003 

00004 

00005 

00006 

00007 

00008 

00009 

00010 
00011 
00012 

00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 
00022 

00023 

00024 

00025 

00026 

00027 

00028 

00029 

00030 

00031 

00032 

00033 

00034 

00035 

00036 

00037 

00038 

00039 

00040 

00041 

00042 

00043 

00044 

00045 

00046 

00047 

00048 

00049 

00050 

00051 

00052 

00053 

00054 

00055 

00056 

00057 


C======== = ==========*===***=********=**=====s===x=»***x**«*s«*«**x=s= 

SUBROUTINE MEC 
C 

C Main Execution Control (MEC) Routine 

C 

C 

PARAMETER (M=5000) 

COMMON /CHECK/ ICHECK.L 

COMMON /ITYPE / ITYPE 

COMMON /PLOTS/ IPLOT 

COMMON /LAST1/ TLAST 

COMMON /FL0W1/ ISW, N, T, DMPDT, T1P0, PIP 

COMMON /FL0W2/ Z1TIME, Z1DMDT, Z1T1P0, Z1P1P, TSTEP 

DIMENSION T (M) , DMPDT (M) , T1P0 (M) , PI P (M) , THRUST (M) , PHI (M) 

C 

C — LOOP ON SPC ROUTINE FOR SPECIFIED PRIMARY FLOW DRIVER 
MODE = 0 
TLAST « T (1) 

DO 20 1=1, N 

IF (I.EQ.l) 60 TO 20 
MODE = 1 
TLAST = T(I-l) 

IF (ISW.EQ.l) GO TO 5 
TSTEP = T(I)-T(I-1) 

GO TO 15 

C — STEP FUNCTION OPERATION — 

5 T(I) = T(I-l) + TSTEP 
IF (T(I).GE.ZITIME) GO TO 6 

T1P0(I) = T1P0(1) 

P1P(I) = P1P(1) 

DMPDT (I) = DMPDT (1) 

GO TO 15 

6 T1P0(I) = Z1T1P0 

P1P(I) = Z1P1P 
DMPDT (I) = Z1DMDT 

GO TO 15 

15 IF (ITYPE. EQ.O) M0DE=0 

20 CALL SPC (T (I ) , DMPDT ( I ) , T1 PO (I) ,P1P(I) , THRUST ( I ) , PHI ( I ) .MODE) 

C 

C — PRINT TO TEXT AND PLOT FILES AS REQUESTED 
IF (ICHECK.GE.50) THEN 
WRITE(6,100) 

WRITE(6, ' ((2X,4(2X,F10.5))) ') 

4 (((T(I)*1000.0) ,DMPDT(I) ,THRUST(I) ,PHI(I)),I=1,N,L) 

END IF 

IF (IPLOT. GT.O) THEN 
WRITE(10, 109) 

WRITEilO, ' ((2X,4(2X,F10.5))) ') 

4 (((T(I)*1000.0), DMPDT(I),THRUST(I),PHI(I)), 1=1, N, IPLOT) 

END IF 
C 

100 FORMAT (1H1/ , 14X,5(1H*) , ' SOLUTION PROFILE ',5(111*), 

4 //7X, ' T ',4X,' DMpDT ' ,6X, 'THRUST* ,8X, 'PHI' , 

4 /7X, ' (m-sec) ' ,4X, ' (lbm/s) ' ,6X, ' (lb)',/) 

109 FORMAT (7X, ' T ' ,4X, ' DMpDT ' ,6X, 'THRUST* ,8X, 'PHI') 

199 RETURN 
END 
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00001 

00002 

00003 

00004 

00005 

00006 

00007 

00008 

00009 

00010 
00011 
00012 

00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 
00022 

00023 

00024 

00025 

00026 

00027 

00028 

00029 

00030 

00031 

00032 

00033 

00034 

00035 

00036 

00037 

00038 

00039 

00040 

00041 

00042 

00043 

00044 

00045 

00046 

00047 

00048 

00049 

00050 

00051 

00052 

00053 

00054 

00055 

00056 

00057 


^xssEsssss9sitigcssisBSStis£>i£s:sssss:sse:s:sts3cs:s:; 

SUBROUTINE SRC (T , DMPDT , T1P0 , PI P t THRUST , PHI .MODE) 
C 

C Single Point Calculation (SPC) Routine 


COMMON /CHECK/ ICHECK, INDEX 
COMMON /COEFX/ BETA, CSIGMA 


COMMON /LAST1/ TLAST 
COMMON /LIST1/ XO, 


COMMON /LIST2/ 6AMMA, 
COMMON /LIST3/ ALPHA, 
COMMON /LIST4/ E, 
COMMON /LISTS/ PIS, 
COMMON /LIST6/ PE, 
COMMON /LIST7/ PM, 
COMMON /LIST8/ DRDT, 


YO, 

Yl, 

Y2, 

Y3, 

Y4, 

NZ1, 

RBAR, 

U 

F. 

TIS, 

TE, 

TM. 

DEDT, 


DZ 

TINF, 

G 

VIS, 

VE, 

VM, 


DIMENSION 

& G(15, 50) , PE (50) , PM (50) , RE(50), RM(50), 

& TE(50), TM(50), VE(50), VM(50), XIHAT(50), Z(50), 

& BI (50) , BI I (50) , DRDT(50) , DPDT(50), DEDT (50) , DVDT(50) 

REAL*8 V1P,V1S,R2,THETA1,DZ,B,BI,BII,THETA2,XIHAT 
REAL*8 SURGE , DELTA, DM3DT , DM4DT , D1S , DIP , A1S , A1PT 
REAL*8 T1S,CIS,PE,PM,VE,VM,RE,RM, DMSDT 
REAL*8 BVALUE,CVALUE,B24AC,C8 
REAL*8 DM1DT, FLUXP, FLUXM, FLUXMI .FLUXMJ 
C 

C — INTERNAL FUNCTION DEFINITIONS 
C 

CP(G,R) - R*G/ (G-1.0) 

F2(G,Z) * 1.0 + 0.5*(G-1.0)*Z*Z 

F3(G,Z) * F2 (G , Z) ** (G/ (G- 1 . 0) ) 

FZ0(PM,PE,E1,E2) * PM*E1 + PE*E2 
FZ1(RM,RE, VM, VE,C1,C2,C3,C4) * 

& RM*VM*C1 + RE*VM*C2 + RM*VE*C3 + RE*VE*C4 
FZ2(RM,RE,VM,VE,G1,G2,G3,G4,G5,G6) = 

& RM*G1*VM*VM + RM*G2*VE*VE + RM*G3*VE*VM + 

& RE*G4*VM*VM + RE*G5*VE*VE + RE*G6*VE*VM 

FZETA(Z,C1,C2,C3,C4) = C1*Z**4. + C2*Z**2.5 + C3*Z + C4 
FPZETA(Z,C1,C2,C3) = 4.*C1*Z**3. + 2.5*C2*Z**1.5 + C3 
FPPZETA(Z,C1,C2) = 12.*C1*Z**2. + 3.75*02*2**0.5 
GFT(T) =1.3825+0. 10173*T*(1.-1.5638*T*(1.0-.4016*T*(1.-.1304*T))) 
H1G(Z) = (4531. + 1309.*Z**10. - 9240.*Z**8.5 + 22440. *Z**7. 

& - 19040. *Z**5. 5)/ 13090.0 

H2G(Z) = (729. + 1122.*Z**10. - 7920.*Z**8.5 + 22440.*Z**7. 

& -32640.*Z**5.5 + 25245. *Z**4. - 8976.*Z**2.5)/ 3740.0 

H3G(Z) = (2727. - 3927.*Z**10. + 27720.*Z**8.5 - 72930. *Z**7. 

& +85680.*Z**5.5 - 39270. *Z**4.)/13090.0 

H4G(Z) = ( 6561. - 2618.*Z**10. + 18480. *Z**8.5 - 56100. *Z**7. 

& +95200.*Z**5.5 - 98175. *Z**4. + 62832.*Z**2.5 

& -26180. *Z)/26180.0 


B(50) , 

PE (50 , 
TM(50), 
811(50), 


PINF, 

PBACK 

D1S, 

Z1S 

RE 


RM 


VENEW, 

, DVE 

50). 

F (15, 50) , 

»). 

RM(50) , 

»). 

XIHAT (50) , Z (50) , 


DEDT (50), DVDT (50) 
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00058 

00059 

00060 
00061 
00062 

00063 

00064 

00065 

00066 

00067 

00068 

00069 

00070 

00071 

00072 

00073 

00074 

00075 

00076 

00077 

00078 

00079 

00080 
00081 
00082 

00083 

00084 

00085 

00086 

00087 

00088 

00089 

00090 

00091 

00092 

00093 

00094 

00095 

00096 

00097 

00098 

00099 

00100 
00101 
00102 

00103 

00104 

00105 

00106 

00107 

00108 

00109 

00110 
00111 
00112 

00113 

00114 


H1L(Z) - -(H1G(Z) - 4531.0/13090.0) 

H2L(Z) = — (H2G (Z) - 729.0/ 3740.0) 

H3L(Z) = — (H3G (Z) - 2727.0/13090.0 
H4L(Z) = -(H4G(Z) - 6561.0/26180.0) 

C 

C — GLOBAL PARAMETER INITIALIZATION 
C 

ITN * 0 
LIMIT = 50 
GC = 32.174 
W = Y2 

SCALE = FL0AT(NX1)*FL0AT(NY1) 

N = NZ1+1 
DZ = Z1/FL0AT(NZ1) 

C 

C — CONDITIONS FOR FREE-STREAM 
C 

10 ALPHAR = ALPHA*3. 1415/180.0 
VINF = U*SIN (ALPHAR) 

DUyDT = DUDT*SIN (ALPHAR) 

DINF = PINF/(RBAR*TINF/GC) 

POINF = PINF+O . 5*V INF*VINF*DINF/GC 
TOINF = TINF+O. 5*VINF*VINF/CP (GAMMA, RBAR) 

AINF = X0*Y0 
DMIDT = AINF*VINF*DINF 
C 

C — EVALUATE PRIMARY NOZZLE CONDITIONS 
C 

A1P = X1*Y1 

A1PT= A1P*FLOAT(NX1)*FLOAT(NY1) 

ZZO = GAMMA/ (GAMMA- 1.0) 

ZZ1 = (GAMMA-1.0) /GAMMA 
ZZ2 = (GAMMA+1.0)/2.0 
ZZ3 = 2 . 0*GAMMA* (GAMMA+1 . 0) 

ZZ4 = (DMPDT/ (A1PT*PBACK*GC)) **2.0 

ZZ4 = ZZ4*T1P0*RBAR* (GAMMA- 1 .0)/(2.0*GAMMA) 

PRN = (0. 5+0. 5*SQRT (1 .0+4.0*ZZ4) ) **ZZO 
XPRN= ZZ2**ZZ0 
C 

C CHECK FOR CHOKED FLOW 

IF (PRN.LT.XPRN) GO TO 17 
C 

C CHOKED NOZZLE FLOW OPERATIONS... 

P1P0 = 1. 88076* (DMPDT/A1PT)*SQRT(T1P0) 

PRN = P1P0/PBACK 

ZIP = 1.0 

PIP = P1P0/F3 (GAMMA, ZIP) 

YPRN = P1P0/P1P 

F0RCE1P= DMPDT* (2 . O/GC) *SQRT ( (T1P0*RBAR)/ZZ3) 
F0RCE1P* F0RCE1P* (GAMMA + 1.0 - XPRN/PRN) 

GO TO 19 


C 


17 


SUBSONIC NOZZLE COMPUTATIONS... 

ZZ4 = (DMPDT/ (A1PT*P1P*GC))**2.0 

ZZ4 = ZZ4*T1P0*RBAR*(GAMMA-1.0)/(2.0*GAMMA) 

YPRN = (0 . 5+0 . 5*SQRT ( 1 . 0+4 . 0*ZZ4) ) **ZZ0 

ZIP = SQRT( (2.0/ (GAMMA-1. 0))*(YPRN**ZZ1 - 1.0)) 
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00115 

00116 

00117 

00118 

00119 

00120 
00121 
00122 

00123 

00124 

00125 

00126 

00127 

00128 

00129 

00130 

00131 

00132 

00133 

00134 

00135 

00136 

00137 

00138 

00139 

00140 

00141 

00142 

00143 

00144 

00145 

00146 

00147 

00148 

00149 

00150 

00151 

00152 

00153 

00154 

00155 

00156 

00157 

00158 

00159 

00160 
00161 
00162 

00163 

00164 

00165 

00166 

00167 

00168 

00169 

00170 

00171 


P1P0 = P1P*YPRN 

PRN = P1P0/PBACK 

PRN = AMAX1(1. ,PRN) 

F0RCE1P® DMPDT* (SQRT (T1P0) ) /GC 

F0RCE1P® FORCE1P*SQRT((2.*RBAR*ZZO)*(1.0-(1./PRN)**ZZ1)) 

19 CONTINUE 

TIP = T1P0/F2 (GAMMA, ZIP) 

DIP = P1P/(T1P*RBAR/GC) 

C1P = SQRT (GAMMA*RBAR*T1P) 

VIP « Z1P*C1P 
DMPDT = D1P*V1P*A1PT 
C 

C — PREDICT SECONDARY FLOW PROPERTIES 
C 

A1S * (X3*Y3)-A1PT 

P1S0*P0INF 

T1S0*T0INF 

C STATION 1-S INITIAL VELOCITY ESTIMATE 
IF (MODE.EQ.O) VIS = 1.0 
IF (MODE.EQ.l) VIS = VENEW 
C ITERATE ON ACTUAL SOLUTION 

20 ITN * ITN + 1 

T1S = T1S0 - 0.5*(V1S**2.0)/CP(GAMMA,RBAR) 

PIS = P1S0*(T1S/T1S0)** (GAMMA/ (GAMMA-1 .0)) 

D1S = PIS/ (RBAR*T1S/GC) 

DMSDT * D1S*V1S*A1S 
ER = DMSDT /DMPDT 
ENERGY1S - 0.5*D1S*V1S**3.0 
C 

C — MIXING REGION PREDICTIONS 
C 

AREA3 = X3*Y3 
DM1DT = DMPDT + DMSDT 
PM(1) = PIP 
PE(1) = PIS 
RM(1) « D1P/GC 
RE(1) - D1S/GC 
VM(l) = VIP 
VE(l) = VIS 

c 

c — VIRTUAL GRID INITIALIZATION 
C 

R2 = V1S/V1P 

THETA1 = DATAN(0.4*(1-R2)*(0.416+0.134*R2)/ (1+R2)) 

BI(1) = A1P/(W*2.0) 

BI(2) >= BI(1) - DZ*DTAN(THETA1) 

IF(BI (2) .LE.O.O) BI(2)=0.0 
BII(l) * BI(1) 

BMAX * X2/2.0 

C 

C COMPUTE OUTER JET EXPANSION ANGLE 

Cl = BI (2)*D1P*V1P + (BMAX-BI (2) ) *D1S*V1S 
C2 = 0.45*D1P*V1P + 0.55*D1P*V1S - D1S*V1S 
B(2) = ( (DM1DT / (2 . 0*SCALE*W) ) -Cl ) /C2 

BII (2) = B(2)+BI(2) 

THETA2= DATAN ( (BII (2)-BII (1) ) /DZ) 

C 
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00172 

00173 

00174 

00175 

00176 

00177 

00178 

00179 

00180 
00181 
00182 

00183 

00184 

00185 

00186 

00187 

00188 

00189 

00190 

00191 

00192 

00193 

00194 

00195 

00196 

00197 

00198 

00199 

00200 
00201 
00202 
00203 

00231 

00232 

00233 

00234 

00235 

00236 

00237 

00238 

00239 

00240 

00241 

00242 

00243 

00244 

00245 

00246 

00247 

00248 

00249 

00250 

00251 

00252 

00253 

00254 

00255 


C COMPUTE FINAL GRID BOUNDARIES 
DO 34 I = 2,N 

Z(I) = DZ*FL0AT(I-1) 

BI(I) * BI(1) - Z(I)*DTAN(THETA1) 

IF(BKI).LE.O.O) BI(I)=0.0 
BII(I) = 81 1(1) + Z(I)*DTAN(THETA2) 

IF(BII (I) .GT.BMAX) BII(I)=BMAX 
8(1) = BII(I) - BI (I) 

XIHAT(I) = (BHAX-BI (I))/B(I) 

34 CONTINUE 
C 

C COMPUTE SIGMA 

C SIGMA = CSIGMA*Z(N)*144.0/(1.0+Z1P) 

SIGMA = CSIGMA*( ((1.0+0. 23*Z1P) *12. 0/Z(N)) **2.0) 

C 

C COMPUTE SELF-SIMILAR PROFILE INTEGRALS FOR THE NEW JET 
DO 37 1=2, N 

BIB = BI(I)/B(I) 

E (1 , 1 ) = BIB + 1.00 
E(I ,2) = XIHAT(I) - 1.00 
F (1 , 1 ) = BIB + 0.45 
F(I,2) = 0.0 
F ( I , 3) = 0.55 
F(I,4) = XIHAT(I) - 1.0 
G(I,1) = (BIB + 243.0/770.0) 

G(I,2) = (320.0/770.0) 

G(I,3) = (414.0/1540.0) 

G(I .4) = 0.0 

G(I,5) = (XIHAT(I) - 1.0) 

G(I,6) = 0.0 
37 CONTINUE 

INCLUDE 'XCHECK1.F0R* 

C 

C SKIP FINITE VOLUME INITIALIZATION IF MODE > 0 

C 

IF (MODE.EQ.l) GO TO 40 
C 

C INITIALIZATION OF FIELD VARIABLES 

C 

DO 39 I = l.N-l 
J = 1+1 
RE (J) = RE(I) 

PE(J) = PE (I) 

DO 38 K-1,3 
VE(J) = VE(I) 

C8 = 0.5*VE(I)*VE(I) 

& + GAMMA/ (GC*(GAMMA-1.0))*(PE(I)/RE(I)-PE(J)/RE(J)) 

IF(C8.GT.0.0) VE(J) = SQRT(2*C8) 

IF (I.EQ.l) THEN 

FLUXRI = RM(I)*VM(I)*A1PT + RE(I)*VE(I)*A1S 
FLUXPI = PM(I)*VM(I)*A1PT + PE(I)*VE(I)*A1S 
FLUXMI = (PM(I)+RM(I)*VM(I)*VM(I))*A1PT 

& + (PE(I)+RE(I)*VE(I)*VE(I))*A1S 

ELSE 

FLUXRI = 2 . 0*SCALE*W*B ( I ) * 

& FZ1(RM(I) ,RE(I) ,VM(I) ,VE(I) ,F(I,1) ,F(I,2) ,F(I,3),F(I,4)) 

FLUXPI = 2. 0*SCALE*W*B ( I ) * 
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00256 

00257 

00258 

00259 

00260 
00261 
00262 

00263 

00264 

00265 

00266 

00267 

00268 

00269 

00270 

00271 

00272 

00273 

00295 

00296 

00297 

00298 

00299 

00300 

00301 

00302 

00303 

00304 

00305 

00306 

00307 

00308 

00309 

00310 

00311 

00312 

00313 

00314 

00315 

00316 

00317 

00318 

00319 

00320 

00321 

00322 

00323 

00324 

00325 

00326 

00327 

00328 

00329 

00330 

00331 

00332 

00333 


& FZKPMdJ.PEdJ.VHdJ.VEdJ.Fd.D.Fd^J.Fd.Sj.Fd^)) 

FLUXMI = 2 . 0*SCALE*W*B ( I ) * ( 

& FZ2(RM(I) ,RE(I) ,VM(I) t VE(I) ,6(1,1) ,G(I,2),6(1,3) ,6(1,4) , 

& G(I,5),6(I,6)) 

& +FZO(PM(I) ,PE(I),E(I,1) |E(I,2))) 

ENDIF 

Cl = (G ( J , 3) +G ( J , 6) ) / (G ( J , 1 )+G ( J , 4) ) 

C2 = (F(J,1)+F(J,2))/(G(J,1)+G(J,4)) 

C3 = (G(J,2)+G(J,5))/(G(J,1)+&(J,4)) 

C4 
C5 


(E(J,1)+E(J,2))/(G(J,1)+G(J,4)) 

.4)) 


(f(j,3)+f(j,4))/(g(j,i)+g(o. 

BVALUE= VE(J)*C1 - (FLUXMI/FLUXRI)*C2 
CVALUE= VE(J)*VE(J)*C3 + (FLUXPI/FLUXRI)*C4 
& - VE(J)*(FLUXHI/FLUXRI)*C5 

B24AC = (BVALUE**2.0) - 4.0*CVALUE 
IF (B24AC.LT. 0.0) THEN 
WRITE d, 900) 1 
INCLUDE •XCHECK2.F0R’ 

ENDIF 

VM(J) * -0.5*(BVALUE+DSQRT (B24AC) ) 

C6 = FLUXRI/(2.0*SCALE*W*B(J)*(F(J,1)+F(J,2))) 

C7 = VE(J)*(F(J,3)+F (J,4) )/(F(J, 1) + F(J«2)) 

C8 = FLUXPI/(2.0*SCALE*W*B(J)*(F(J,1)+F(J,2))) 
RM(J) = C6/(VM(J)+C7) 

RE(J) = RM(J) 

PM(J) - C8/(VM(J)+C7) 

PE(J) = PM(J) 

38 CONTINUE 

39 CONTINUE 


C — DERIVATIVE ESTIMATES 
C 

40 DT = T-TLAST 


C 

C ...SKIP DERIVATIVE COMPUTATIONS IF DT-0 
IF(DT.LE.O.O) GO TO 43 
C 

C ...NO DERIVATIVES ARE NEEDED FOR QUASI-STEADY FLOW 
IF(MODE.EQ.O) GO TO 43 
C 

C ...PROCEED WITH DERIVATIVE COMPUTATIONS 
DO 42 I = 1,N-1 
J = 1+1 
C 

C ...UPDATE ENTRAINED VELOCITY 

VE(J) - VE(I) 

C8 = 0.5*VE(I)*VE(I) 

& + GAMMA/ (GC* (GAMMA-1. 0))*(PE(I)/RE(I)-PE(J)/RE(J)) 

IF(C8.GT.0.0) VE(O) = SQRT(2.0*C8) 


C 

C ...COMPUTE STATION 'I' FLUXES 

IF (I.EQ.l) THEN 

FLUXRI = DM1DT/(2.0*SCALE*W*GC) 

FLUXMI = ( ( A1PT* ( PI P+Dl P*V 1 P*V 1 P/GC) ) 

& +(A1S*(P1S+D1S*V1S*V1S/GC)))/(2.0*SCALE*W) 

FLUXPI = ((A1PT*P1P*V1P)+(A1S*P1S*V1S))/(2.0*SCALE*W) 
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00334 

00335 

00336 

00337 

00338 

00339 

00340 

00341 

00342 

00343 

00344 

00345 

00346 

00347 

00348 

00349 

00350 

00351 

00352 

00353 

00354 

00355 

00356 

00357 

00358 

00359 

00360 

00361 

00362 

00363 

00364 

00365 

00366 

00367 

00368 

00369 

00370 

00371 

00372 

00373 

00374 

00375 

00376 

00377 

00378 

00379 

00380 

00381 

00382 

00383 

00384 

00385 

00386 

00387 

00388 

00389 

00390 


C 

C 


C 

C 


C 

C 


C 

C 


ELSE 

FLUXRI = B(I)*FZ1(RM(I) ,RE(I) , VM(I) ,VE(I), 

& F(I,l),F(I,2),F(l,3),F(I,4)) 

FLUXMI = B(I)*(FZ2(RM(I) ,RE(I) ,VH(I) ,VE(I) , 

& G(I,1),6(I|2) t G(I,3) l G(I,4), 

& 6(1.5) ,6(1,6)) 

& +FZO(PM(I) ,PE(I j , E(I,1),E(I,2))) 

FLUXPI = B(I)*FZl(PM(I),PE(I),VM(I),VE(I) f 
4 F(I, 1) ,F(I,2) ,F(I,3) ,F(I,4)) 

ENOIF 

... COMPUTE STATION * J * FLUXES ... 

FLUXRO = B(J)*FZ1(RM(J),RE(J),VM(J),VE(J), 

& F(J,l),F(J l 2),F(J,3),F(J,4j) 

FLUXMJ = B(J)*(FZ2(RM(J) ,RE(J) , VM(J) ,VE(J) , 

& 6(J, 1) ,G(J,2) ,G(J,3) ,G(J,4) , 

& G(J,5),G(J,6)) 

& +FZO(PM(J) ,PE(J) ,E(J,1) ,E(J,2))) 

FLUXPJ = B(J)*FZ1(PM(J),PE(J),VM(J),VE(J), 

& F(J,1) , F(J,2) ,F(J,3) ,F(J,4)) 

... COMPUTE PRIMARY FLOW DERIVATIVES ... 

DEDT(I) = DVE/DT 

DRDT(I) = (FLUXRI - FLUXRJ)/(BMAX*DZ) 

DPDT(I) = (FLUXPI - FLUXPJ) /(BMAX*DZ) 

CDVDT =0.0 

Cl = (F(J,3)+F(J,4))/ (F(J,1)+F(J,2)) 

C2 = C1*DEDT ( I ) *CDVDT - (DRDT(I)/RM(J))*(VM(J)+VE(J)*C1) 
DVDT(I) = (FLUXMI - FLUXMJ)/(B(J)*DZ) - C2 

... ADVANCE PRIMARY FLOW FORWARD IN TIME ... 

RM(J) = RM(J) + DRDT(I)*DT 
RE(J) = RM(J) 

PM(J) = PM(J) + DPDT(I)*DT 
PE(J) = PM(J) 

VM ( J ) = VM(J) + DVDT(I)*DT 

42 CONTINUE 

COMPUTE THE JET STREAMLINE 

43 DO 45 J = 2,N 

ISTOP= 0 

ZETA = BI(1)/B(J) 

Cl = 0. 25* (VM(J) -VE(J) ) 

C2 = 0.80*(VE(J)-VM(J) ) 

C3 = VM(J) 

C4 = (RE(1)/RE(J))*((BMAX-BI(1))/B(J))*V1S 
& - Cl - C2 - C3 - ((BMAX/B(J)) - l.)*VE(J) 

44 ISTOP = ISTOP +1 

C5 = FZETA(ZETA t Cl,C2,C3,C4)/FPZETA(ZETA,Cl,C2,C3) 

C6 = 1- FZETA(ZETA,Cl,C2,C3,C4)*FPPZETA(ZETA t Cl,C2) 

& /FPZETA(ZETA,C1,C2,C3)**2.0 

DELTA = -C5/C6 
ZETA = ZETA + DELTA 
BSNEW = ZETA*B(J) 

IF (BSNEW. GE.BMAX) THEN 
BSTAR(J) = BMAX 
GO TO 45 
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00391 

00392 

00393 

00394 

00395 

00396 

00397 

00398 

00399 

00400 

00401 

00402 

00403 

00404 

00405 

00406 

00407 

00408 

00409 

00410 

00411 

00412 

00413 

00414 

00415 

00416 

00417 

00438 

00439 

00440 

00441 

00442 

00443 

00444 

00445 

00446 

00447 

00448 

00449 

00450 

00451 

00452 

00453 

00454 

00455 

00456 

00457 

00458 

00459 

00460 

00461 

00462 

00463 

00464 

00465 

00466 

00467 


ENDIF 

IF ( (ABS (DELTA) .GE. 0.000001). AND. (ISTOP.LT. 50)) GO TO 44 
BSTAR(J) « BSNEW 
45 CONTINUE 
C 

C — KINETIC ENERGY EXCHANGE 
C 

ZETA = BSTAR(N)/B(N) 

ZSTAR - BMAX/B(N) 

Cl = (HIG(ZETA) - (0.55 + 0.25*ZETA**4.-0.8*ZETA**2.5)) 

& *VE(N)**3.0 

C2 = (H3GIZETA) - (0.45-(0.25*ZETA**4.-0.8*ZETA**2.5+ZETA))) 
& *VM(N)*VE(N)**2.0 

C3 « H2G(ZETA)*VE(N)*VM(N)**2.0 
C4 « H4G(ZETA)*VM(N)**3.0 
DKE = B(N)*RM(N)*GC*(C1+C2+C3+C4) 

Cl = H1L(ZETA)*VE(N)**3. + H2L(ZETA)*VE(N)*VM(N)**2. 

& + H3L(ZETA)*VM(N)*VE(N)*2. + H4L(ZETA)*VH(N)**3. 

C2 = (0.25*ZETA**4. - 0.8*ZETA**2.5 + ZETA)*VH(N)**3. 

& + (-0.25*ZETA**4. + 0.8*ZETA**2.5)*VE(N)*VM(N)**2. 

PLOSS * -B(N)*RM(N)*GC*(C1-C2) 

DKE = DKE + SIGMA* PLOSS 

CHANGE - DKE - DKEOLD 

VSTAR = ( (ENERGY1S+CHANGE)/ (0.5*D1S)) 

IF(VSTAR.LT.O.O) THEN 
WRITE (*,900) 2 
INCLUDE 'XCHECK3.F0R 1 
ENDIF 

VENEW * (ABS (VSTAR) )**(l./3.) 

DVE = VENEW-V1S 
DKEOLD = DKE 
C 

C — ASSIGN MIXING REGION EXIT CONDITIONS 
C 

ZZZ = (2 . 0*SCALE*W*B (N) ) /AREA3 

V3 = ZZZ*(0.45*VM(N) + (XIHAT(N)-0.45)*VE(N)) 

D3 « RM(N)*GC 

IF (D3.LT.0.0) WRITE (*,900) 3 
IF (D3.GE.1.0) WRITE (*,900) 4 
P3 = PM(N) 

T3 = P3/(RBAR*D3/GC) 

IF (T3.LT.0.0) WRITE (*,900) 5 
C3 = SQRT (GAMMA*RBAR*T3) 

Z3 = V3/C3 

P30 = P3*F3 (GAMMA, Z3) 

T30 = T3*F2 (GAMMA, Z3) 

DM3DT= D3*V3*AREA3 
C 

C — DIFFUSER PREDICTIONS 
C 

50 AREA4 * X4*Y4 
P40 = P30 

T40 = T30 

IF(MODE.EQ.O) THEN 
P4 = PBACK 

IF(P40.LT.P4) WRITE(*,900) 6 
Z4 = SQRT ((2.0/ (GAMMA-1.0)) 
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00481 

00482 
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00484 

00485 

00486 

00487 

00488 

00489 

00490 

00491 

00492 

00493 

00494 

00505 

00506 

00507 

00508 

00509 

00510 

00511 

00512 

00513 

00514 

00515 

00516 

00517 

00518 

00519 

00520 

00521 
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00524 

00525 

00526 

00527 

00528 
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00530 
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& * ( ( (P40/P4) ** ( (GAMMA-1 . 0) /GAMMA) ) -1 .0) ) 

T4 = T40/F2 (GAMMA, Z4) 

04 = P4/ (RBAR*T4/GC) 

C4 = SQRT (GAMMA*RBAR*T4) 

V4 = Z4*C4 
ELSE 

Cl = 0.5/ (GAMMA-1.0) 

C2 = (GAMMA/ (GAMMA-1. 0))*GC*AREA4*P30/DM3DT 

C3 = CP(GAMMA,RBAR)*T30 

V4 * (C2-SQRT(C2*C2-4.0*C1*C3))/(2.0*C1) 

04 = DM3DT/(V4*AREA4) 

P4 = P40 - 0.5*D4*V4*V4/GC 

T4 = T40 - 0.5*V4*V4/CP(GAMMA,RBAR) 

C4 = SQRT (GAMMA*RBAR*T4) 

ENDIF 

DM4DT = D4*V4*AREA4 
C 

C — CHECK FOR CONTINUITY 
C 

EPSILON = 0.00005 

SURGE = (DM4DT - DM3DT)/(D1S*A1S) 

IF(MODE.EQ.l) GO TO 70 
IF (ABS (SURGE). LE. EPSILON) THEN 
GO TO 70 
ELSE 

VIS = VIS + SURGE 
INCLUDE 'XCHECK4.F0R' 

IF(ITN.LT. LIMIT) GO TO 20 
ENDIF 

C THRUST AND THRUST AUGMENTATION RESULTS 


70 THRUST = (DM4DT* ( V4/ (BETA**2. ) -VINF) / (GC) ) 
& + AREA4* (PBACK-PINF) 

PHI = THRUST/F0RCE1P 
C 

C — PRINT OPTION 
C 


IF ( ICHECK.LT. 100) GO 

TO 999 



WRITE (6,200) 

T 




WRITE (6,201) 

AINF, 

A1PT, 

A1S, 

AREA3, AREA4 

& 

VINF, 

VIP, 

VIS, 

V3, V4, 

& 

PINF, 

PIP, 

PIS, 

P3 ■ P4, 

P30, P40, 

& 

POINF, 

P1P0, 

P1S0, 

& 

TINF, 

TIP. 

T1S, 

T3, T4, 

& 

TOINF, 

T1P0, 

T1S0, 

T30, T40, 

& 

DINF, 

DIP, 

D1S, 

D3, D4, 

& 

OMIDT, 

DMPDT, 

DMSDT, 

DM3DT, DM4DT 

WRITE (6,202) 

PRN, 

XPRN 




IF(MODE.EQ.O) 
WRITE (6,204) 


WRITE (6,203) SURGE, ITN, LIMIT 


WRITE (6,207) ( ( I ,Z ( I) , BI (I ) , BI I ( I) , BSTAR(I) ) , I=2,N) 

WRITE (6,206) 

WRITE (6,207) ((I,RM(I) ,PM(I) , VM(I) , VE(I)) , 1=1, N) 

WRITE (6,208) 

WRITE (6,207)((I,0RDT(I),DPDT(I),DVDT(I),DEDT(I)),I=1,N-1) 
WRITE (6,209) THETA1,THETA2,BI(1) ,BMAX 
WRITE (6,210) THRUST, F0RCE1P, PHI 
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APPENDIX B - TEA PROGAM LISTING 


00535 200 FORMAT (1H1./.19X, ****** SIN 6 LE POINT CALCULATION *****’, 

00536 4 /,29X, 'T = ',F8.5,' sec'/) 

00537 201 FORMAT ( /, 3 X, 'FIELD VARIABLE PROFILE:',/ 

00538 & /,3X, 'Station -> Infinity', 6 X, 'IP', 10X, 'IS', 10X, 

00539 & 1 3 4 1 / 

00540 & /,3X, 'AREA ' ,2X,5(F10.5,2X) , 'FT' , 

00541 & /,3X, ' V ' , 2 X, 5(F10.5,2X) , 'FT/S' , 

00542 4 /,3X,'P ' ,2X,5(F10.5,2X) , 'LBF/FT2' , 

00543 & /,3X,'P0 ' ,2X,5(F10.5,2X) , 'LBF/FT2' , 

00544 4/,3X,'T ',2X,5(F10.5,2X), 'DEG-R', 

00545 & /,3X, 'TO ' ,2X,5(F10.5,2X) , 'DEG-R' , 

00546 4/,3X,'RH0 ' ,2X,5(F10.5,2X) , 'LBM/FT3' . 

00547 & /|3X, 'MOOT ' ,2X,5(F10.5,2X) , ' LBM/S ' ) 

00548 202 F0RMAT(7X, 'NOTES: ' ,/10X, 'Nozzle Pressure Ratio Is \F8.5, 

00549 & ' (Choked NPR Is \F8.5, ' )') 

00550 203 FORMAT ( 10X, 'Station 3-4 SURGE error is '.F12.8, 

00551 4/, 10X, 'Solution obtained In ',13, ' Iterations, limit Is *,I3) 

00552 204 F0RMAT(/,3X, 'FINITE VOLUME PROFILE', //,5X, 'Station Z', 

00553 410X, 'B, Inner', 7X,'B, outer', BX,'B,jsl ') 

00554 206 FORMAT ( / , 5 X , ' Stat 1 on ' , 4X , ' Rho ’ , 1 IX , ’ P ' , 13X , ' Vm ’ , 12X , ' Ve ' ) 

00555 207 F0RMAT(6X, I3,2X,E12.5,2X,E12.5,2X,E12.5,2X,E12.5) 

00556 208 F0RMAT(//,5X, 

00557 4 'Element 0 (Rho)/DT D(P)/DT D(Vm)/DT 0(Ve)/DT') 

00558 209 FORMAT (7X , ' NOTES : ' , 

00559 4 /10X, 'Inner jet expansion, Thetal ',F8.5, ' radians', 

00560 4 / 10 X, 'Outer jet expansion, ThetaZ ',F8.5, ' radians', 

00561 4 /10X, 'Initial jet with, BO, \F8.5,' Ft', 

00562 4 /10X, 'Allowable jet width, Bmax, ',F8.5,' Ft') 

00563 210 FORMAT ( /,3X, 'PREDICTED THRUST AND PERFORMANCE:', 

00564 4 // , 5X, 'Total ejector thrust = ',F12,5, ' Lb', 

00565 4 /,5X, 'Isentroplc primary thrust = ',F12.5,' Lb', 

00566 4 /,5X, 'Thrust augmentation ratio * ',F12.5) 

00567 900 FORMAT (/5X,'$$$ TEA FATAL ERROR # ’,12,', $$$'/) 

00568 999 RETURN 

00569 END 
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APPENDIX B - TEA PROGAM LISTING 


00001 

00002 

00003 

00004 

00005 

00006 

00007 

00008 

00009 

00010 
00011 
00012 

00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 
00022 

00023 

00024 

00025 

00026 
00027 


C XCHECK1.F0R 


C 

C- 


PRINT MIXING REGION DATA FOR DEBUG CHECKS 


1000 IF(ICHECK.LT.IOOO) GO TO 1020 

WRITE(6,1010) T, ITN 
WRITE(6,1011) 

DO 1001 I=2,N 

1001 WRITE(6, 1012) (I, (E(I,J) , J=l,2)) 
WRITE(6,1013) 

DO 1002 1=2, N 

1002 WRITE(6, 1014) (I, (F(I,J),J=1,4)) 
WRITE(6,1015) 

DO 1003 1=2, N 

1003 WRITE(6, 1016) (I, (G ( I , J) , J=1 , 6) ) 

1010 FORMAT(1H1,/10X, '*** T = \F10.5,’ ITN « \I5) 

1011 FORMAT (/,33X,'> SELF-SIMILAR PROFILE EVALUATION 

& /,3X, 'Z-POINT E ( 1) E{2) ') 

1012 F0RMAT(5X, I2,1X,2(3X,F7.3)) 

1013 F0RMAT(/,3X, ’Z-POINT F(l) F(2) F(3) 

& * F (4) ' ) 

1014 FORMAT (5X, I2,1X,4(3X, F7.3)) 

1015 F0RMAT(/,3X, 'Z-POINT G(l) G(2) G(3) 

& * G (4) G (5) G(6) 

1016 FORMAT (5X, I2,1X,6(3X, F7.3)) 

1020 CONTINUE 


SSSSXSSSSSXESSSSS 


<\ 


•) 


I 

f 
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APPENDIX B - TEA PROGAM LISTING 


00001 

00002 

00003 

00004 

00005 

00006 

00007 

00008 

00009 

00010 
00011 
00012 

00013 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 


C XCHECK2.F0R 

C - 

C PRINT VELOCITY INITIALIZATION DATA FOR DEBUG CHECKS 

C 

2000 IF(ICHECK.LT. 500) GO TO 2060 
WRITE(6,2010) T, I, J 
WRITE (6, 2020) C1,C2,C3,C4,C5 
WRITE (6 ,2050) THETA1 t THETA2 
WRITE (6 , 2030) FLUXRI , FLUXHI , FLUXPI 
WRITE(6,2040) VE(J) ,BVALUE,CVALUE,B24AC 
2010 F0RMAT(///,5X, '*** T = ' ,F7.4, ' ; I = ',12,' J = *,12) 

2020 FORMAT(/ t 5X, 'Cl' ,8X, 1 C2 1 ,8X, 'C3' ,8X, 'CA'.BX, ’C5' , 
&/,5X,5(F10.4,2X)) 

2030 FORMAT (/2X, 1 FLUXRI = \E12.4, 

& /2X, 'FLUXMI = \E12.4,/2X,’ FLUXPI = ’,£12.4) 

2040 F0RMAT(/2X, 'VE(J) = '.F10.4,/ 

& 2X,*B « ' ,E12.4, ' C * \E12.4,' B24AC * \E12.4) 

2050 FORMAT (/5X,’THETA1 = \F10.5,’ THETA2 = \F10.5) 

2060 CONTINUE 
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APPENDIX B - TEA PROGAM LISTING 

00002 C XCHECK3.F0R 

00003 C 

00004 C KINETIC ENERGY EXCHANGE ERRORS 

00005 C Print out all kinetic energy Information If the program bombs. 

00006 C This error by resetting the KE-CHANGE and re-computing VSTAR 

00007 c : 

00008 WRITE (*,*) 

00009 WRITE(*, ' (/5X.38H*** ERROR DURING CALCULATIONS FOR T = , 

00010 & F10.5) ') T 

00011 WRITE (V(/10X,10HDKE = ,F15.3,/10X,10HDKEOLD * .F15.3, 

00012 & /10X.10HSIGMA = ,F15.3,/10X,10HPL0SS = ,F15.3, 

00013 & /10X t 10HENERGY1S= ,F15.3,/10X,10HCHANGE = ,F15.3, 

00014 & /10X, 10HVSTAR = ,F15.3)') 

00015 & DKE.DKEOLD, SIGMA, PL0SS.ENERGY1S, CHANGE, VSTAR 

00016 CHANGE = -0.999*ENERGY1S 

00017 WRITE (*, '(/5X.38H*** PROCEED W/ CHANGE = -0.999ENERGY1S) *) 

00018 VSTAR - ( (ENERGY1$+CHANGE)/(0.5*D1S)) 

00019 WRITE(*,'(/5X,19H*** NEW VSTAR IS « .F15.3//)') VSTAR 
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APPENDIX B - TEA PROGAM LISTING 

00001 ^:s:s»««ii«iK»xs::sss::::ssr:essissXBs:ss:s;ssss«Etistiss3s:: 

00002 C XCHECK4.F0R 

00003 C 

00004 C PRINT ENTRAINED VELOCITY AND SURGE FOR ITERATIVE PROCESS 

00005 C 

00006 4000 IF(ICHECK.LT.IOOO) GO TO 4020 

00007 WRITE(6,4010) VIS, SURGE 

00008 4010 FORMAT (/,5X, 'VIS = * , F10 . 3 f * SURGE = 1 ,F10.3) 

00009 4020 CONTINUE 

00010 C 


32 



APPENDIX C - TEST CASE INPUT LISTINGS 


C - CASE1.DAT — 

C 

C Case study Input file for TEA. FOR transient ejector analysis 

C routine. Here, a single step function is prescribed that 

C corresponds with one of the STOVL ejector tests run in June 1987. 

C This file assumes the default BACK and DISP parameters In the TEA. FOR 
C source code are appropriate. 

C 

C 


CASE banner. 

Test for step change In primary flow for STOVL ejector design. 

6E0Metr1c approximation for the ejector. 

XO YO 

0009.36000, 0001.25000 

XI Y1 NX1 NY1 

0000.10801, 0000.10801, 00012, 00003 

X2 Y2 

0000.78000, 0000.34700 

X3 Y3 

0009.36000, 0001.04100 

X4 Y4 

0009.36000, 0001.87500 

Z1 NZ 

0000.90000, 00005 

COEFfl clents for solution execution. 

ITYPE BETA CSIGMA 

00001, 0001.0550, 0000.30000 

PRIMary nozzle discharge conditions. 

SPECIfy states at each time step. 

T MDOT1P T1P0 PIP 

0000.00010, 0018.70000, 0769.70000, 2219.78000 
0000.00020, 0018.70000, 0769.70000, 2219.78000 
0000.00030, 0018.70000, 0769.70000, 2219.78000 
0000.00040, 0018.70000, 0769.70000, 2219.78000 
0000.00050, 0018.70000, 0769.70000, 2219.78000 
0000.00060, 0021.85000, 0760.60000, 2250.92000 
0000.00070, 0021.85000, 0760.60000, 2250.92000 
0000.00080, 0021.85000, 0760.60000, 2250.92000 
0 

C 
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00001 

00002 CASE2.DAT — 

00003 C , J 

00004 C This case study has the same basic conditions as CASE1.DAT, but Includes 

00005 C the use of some optional macros. 

00006 C 

00007 C - 

00008 


00009 CASE banner. 

00010 Test for step change In primary flow for STVOL ejector design. 

00011 

00012 DISPlay parameters for printing and plotting output. 

00013 ICHECK INDEX IPLOT 

00014 00100, 00001, 00001 

00015 

00016 6E0Metr1c approximation for the ejector. 

00017 XO YO 

00018 0009.36000, 0001.25000 

00019 XI Y1 NX1 NY1 

00020 0000.10801, 0000.10801, 00012, 00003 

00021 X2 Y2 

00022 0000.78000, 0000.34700 

00023 X3 Y3 

00024 0009.36000, 0001.04100 

00025 X4 Y4 

00026 0009.36000, 0001.87500 

00027 Z1 NZ 

00028 0000.90000, 00005 

00029 

00030 COEFfl clents for solution execution. 

00031 ITYPE BETA CSIGMA 

00032 00001, 0001.0550, 0000.30000 

00033 

00034 FREE stream fluid conditions. 

00035 TINF PINF 6AMMA RBAR UINF ALPHA 

00036 0542.80000, 2067.84000, 0001.40000, 1714.54130, 0000.00000, 0090.00000 

00037 

00038 BACKpressure for ejector discharge 

00039 2067.84000 

00040 

00041 PRIMary nozzle discharge conditions. 

00042 SPECIfy states at each time step. 

00043 T HD0T1P T1P0 PIP 

00044 0000.00010, 0018.70000, 0769.70000, 2219.78000 

00045 0000.00020, 0018.70000, 0769.70000, 2219.78000 

00046 0000.00030, 0018.70000, 0769.70000, 2219.78000 

00047 0000.00040, 0018.70000, 0769.70000, 2219.78000 

00048 0000.00050, 0018.70000, 0769.70000, 2219.78000 

00049 0000.00060, 0021.85000, 0760.60000, 2250.92000 

00050 0000.00070, 0021.85000, 0760.60000, 2250.92000 

00051 0000.00080, 0021.85000, 0760.60000, 2250.92000 

00052 0000.00090, 0021.85000, 0760.60000, 2250.92000 

00053 0000.00100, 0021.85000, 0760.60000, 2250.92000 

00054 0000.00110, 0021.85000, 0760.60000, 2250.92000 

00055 0000.00120, 0021.85000, 0760.60000, 2250.92000 

00056 0000.00130, 0021.85000, 0760.60000, 2250.92000 

00057 QUIT primary nozzle data Input 

00058 

00059 STOP reading DATA 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 


***** tea INPUT ***** 


Test for step change In primary flow for STVOL ejector design. 
CHECK PRINT LEVEL IS 100 


EJECTOR GEOMETRIC CHARACTERISTICS, 


XO = 

9.36000 

YO = 

1.25000 

XI = 

0.10801 

Y1 * 

0.10801 

NX1= 

12 

NY1= 

3 

X2 * 

0.78000 

Y2 = 

0.34700 

X3 = 

9.36000 

Y3 = 

1.04100 

X4 = 

9.36000 

Y4 = 

1.87500 

Z1 = 

0.90000 

NZ = 

5 


MIXING LOSS CORRECTION FACTOR 
BETA = 1.05500 

KINETIC ENERGY MIXING COEFFICIENT 
CSIGMA = 0.30000 

FREE-STREAM THERMODYNAMIC PROPERTIES 
STATIC TEMPERATURE = 542.79999 DEG-RANKINE 
STATIC PRESSURE = 2067.84009 LB-F/FT3 
SPECIFIC HEAT RATIO = 1.40000 

GAS CONSTANT, AIR = 1714.54126 FT2/S2-R 
FREESTREAM VELOCITY = 0.00000 FT/S 

ANGLE-OF-ATTACK = 90.00000 DEGREES 


EJE BACK PRESSURE IS 2067.84009 LBF/FT2 


PRIMARY 

NOZZLE DATA 




DATA 

T 

DMPDT 

T1P0 

PIP 

SET# 

(sec) 

(Ibm/s) 

(degR) 

(lbf/ft2) 

1 

0.00040 

18.70000 

769.70001 

2219.78003 

2 

0.00050 

18.70000 

769.70001 

2219.78003 

3 

0.00060 

21.85000 

760.59998 

2250.91992 

4 

0.00070 

21.85000 

760.59998 

2250.91992 

5 

0.00080 

21.85000 

760.59998 

2250.91992 

PROGRAM 

WILL OPERATE 

ON 5 

DATA SETS 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 

***** SINGLE POINT CALCULATION ***** 

T = 0.00040 sec 


FIELD VARIABLE PROFILE: 





Station 

-> Infinity 

IP 

IS 

3 

4 


AREA 

11.70000 

0.41998 

9.32378 

9.74376 

17.55000 

FT 

V 

0.00000 

769.94739 

282.19669 

304.57248 

164.99725 

FT/S 

P 

2067.84009 

2219.78003 

1980.71228 

1997.96570 

2067.84009 

LBF/FT2 

PO 

2067.84009 

2799.79346 

2067.84009 

2097.37158 

2097.37158 

LBF/FT2 

T 

542.79999 

720.30579 

536.16473 

553.28766 

558.74860 

DEG-R 

TO 

542.79999 

769.70001 

542.79999 

561.01691 

561.01691 

DEG-R 

RHO 

0.07149 

0.05783 

0.06932 

0.06776 

0.06945 

LBM/FT3 

MOOT 

0.00000 

18.70000 

182.39970 

201.09969 

201.09970 

LBM/S 

NOTES: 




1.89293 ) 



Nozzle Pressure Ratio Is 

1.35397 (Choked NPR Is 



Station 3-4 SURGE error Is 0.00002361 
Solution obtained In 17 Iterations, limit Is 50 


FINITE VOLUME PROFILE 


Station Z 

B, Inner 

B, outer 

B,jsl 

2 

0. 18000 E+00 

0.12857E-02 

0.42297E-01 

0.10382E-01 

3 

0.36000E+00 

O.OOOOOE+OO 

0.67784E-01 

0.16039E-01 

4 

0.54000 E+00 

0.00000E+00 

0.93272E-01 

0.19376E-01 

5 

0.72000E+00 

0.00000E+00 

0.11876E+00 

0.21867E-01 

6 

0.90000 E+00 

0.00000E+00 

0.14425E+00 

0.23816E-01 

Station Rho 

P 

Vm 

Ve 

1 

0.17974E-02 

0.22198E+04 

0.76995E+03 

0.28220E+03 

2 

0.20635E-02 

0.19575E+04 

0.10653E+04 

0.27064E+03 

3 

0.20895E-02 

0 . 19822E+04 

0.73541E+03 

0.27064E+03 

4 

0. 20987 E-02 

0.19909E+04 

0.59600E+03 

0.27064E+03 

5 

0.21033E-02 

0.19953E+04 

0.52124E+03 

0.27064E+03 

6 

0.21061E-02 

0.19980E+04 

0.47450E+03 

0.27064E+03 

Element D(Rho)/DT 

D(P)/DT 

D(Vm)/DT 

D(Ve)/DT 

1 

0.00000 E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

2 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

3 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

4 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

5 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 


NOTES: 

Inner jet expansion, Thetal 0.08603 radians 
Outer jet expansion, Theta2 0.14066 radians 
Initial Jet with, B0, 0.01681 Ft 

Allowable jet width, 8max, 0.39000 Ft 

PREDICTED THRUST AND PERFORMANCE: 

Total ejector thrust = 926.57001 Lb 

Isentroplc primary thrust = 508.74969 Lb 

Thrust augmentation ratio * 1.82127 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 

***** SINGLE POINT CALCULATION ***** 

T = 0.00050 sec 


FIELD VARIABLE PROFILE: 

Station -> Infinity IP IS 3 4 

AREA 11.70000 0.41998 9.32378 9.74376 17.55000 FT 

V 0.00000 769.94739 282.19650 304.57236 164.98514 FT/S 

P 2067.84009 2219.78003 1980.71240 1997.96570 2067.99219 LBF/FT2 

PO 2067.84009 2799.79346 2067.84009 2097.37158 2097.37158 LBF/FT2 

T 542.79999 720.30579 536.16474 553.28766 558.74890 DEG-R 

TO 542.79999 769.70001 542.79999 561.01691 561.01691 DEG-R 

RHO 0.07149 0.05783 0.06932 0.06776 0.06945 LBM/FT3 

MOOT 0.00000 18.70000 182.39959 201.09961 201.09961 LBM/S 

NOTES: 

Nozzle Pressure Ratio Is 1.35397 (Choked NPR Is 1.89293 ) 

FINITE VOLUME PROFILE 


Station Z 

B, Inner 

B, outer 

B.jsl 

2 

0.18000E+00 

0.12857E-02 

0.42297E-01 

0.10382E-01 

3 

0.36000E+00 

O.OOOOOE+OO 

0.67784E-01 

0.16039E-01 

4 

0.54000E+00 

0.00000E+00 

0.93272E-01 

0.19376E-01 

5 

0.72000E+00 

O.OOOOOE+OO 

0.11876E+00 

0.21867E-01 

6 

0. 90000 E+00 

O.OOOOOE+OO 

0.14425E+00 

0.23816E-01 

Station Rho 

P 

Vm 

Ve 

1 

0.17974E-02 

0.22198E+04 

0.76995E+03 

0.28220 E+03 

2 

0.20635E-02 

0.19575E+04 

0.10653E+04 

0.27064E+03 

3 

0.20895E-02 

0.19822E+04 

0.73541E+03 

0.27064E+03 

4 

0.20987E-02 

0.19909E+04 

0. 59600 E+03 

0.27064E+03 

5 

0.21033E-02 

0.19953E+04 

0.52124E+03 

0.27064E+03 

6 

0. 21061 E-02 

0.19980E+04 

0.47450E+03 

0.27064E+03 

Element D(Rho)/DT 

D(P)/DT 

D(Vm)/DT 

D(Ve)/DT 

1 

O.OOOOOE+OO 

0.44516E+00 

0.71454E-02 

-0.18450E+01 

2 

0.00000E+00 

0.22258E+00 

0.61080E-02 

-0.18450E+01 

3 

0. 42453 E-06 

0.44516E+00 

0.57823 E+00 

-0. 18450E+01 

4 

-0. 42453 E-06 

0.22258E+00 

-0.444 30 E+OO 

-0.18450E+01 

5 

-0.42453 E-06 

-0.22258E+00 

-0.37205E+00 

-0.18450E+01 


NOTES: 

Inner jet expansion, Thetal 0.08603 radians 
Outer jet expansion, Theta2 0.14066 radians 
Initial jet with, BO, 0.01681 Ft 
Allowable jet width, Bmax, 0.39000 Ft 

PREDICTED THRUST AND PERFORMANCE: 

Total ejector thrust = 926.50159 Lb 

Isentroplc primary thrust = 508.74969 Lb 

Thrust augmentation ratio = 1.82113 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 


***** SINGLE POINT CALCULATION ***** 
T * 0.00060 sec 


FIELD VARIABLE PROFILE: 


Station 

-> Infinity 

IP 

IS 

3 

4 


AREA 

11.70000 

0.41998 

9.32378 

9.74376 

17.55000 

FT 

V 

0.00000 

860.78772 

282.19669 

303.30927 

164.33514 

FT/S 

P 

2067.84009 

2250.91992 

1980.71228 

1998.06897 

2067.49927 

LBF/FT2 

PO 

2067.84009 

3027.15869 

2067.84009 

2096.64307 

2096.64307 

LBF/FT2 

T 

542.79999 

698.86292 

536.16473 

553.28772 

558.70282 

DEG-R 

TO 

542.79999 

760.59998 

542.79999 

560.95300 

560.95300 

DEG-R 

RHO 

0.07149 

0.06044 

0.06932 

0.06777 

0.06944 

LBM/FT3 

MOOT 

0.00000 

21.85001 

182.39970 

200.27596 

200.27597 

LBM/S 


E5: 

Nozzle Pressure Ratio Is 1.46392 (Choked NPR Is 1.89293 ) 


FINITE VOLUME PROFILE 


Station Z 

B, Inner 

B, outer 

B.jsl 

2 

0.18000E+00 

0.46881E-04 

0.41181E-01 

0.11199E-01 

3 

0 . 36000E+00 

0.00000E+00 

0.65551E-01 

0.15560E-01 

4 

0.54000E+00 

0.00000E+00 

0.89922E-01 

0.18624E-01 

5 

0.72000E+00 

0 . 00000E+00 

0.11429E+00 

0.20934E-01 

6 

0.90000E+00 

O.OOOOOE+OO 

0.13866E+00 

0.22756E-01 

Station Rho 

P 

Vm 

Ve 

1 

0. 18785E-02 

0.22509E+04 

0.86079E+03 

0.28220E+03 

2 

0.20719E-02 

0.19656E+04 

0.10916E+04 

0.27064E+03 

3 

0.20910E-02 

0.19836E+04 

0.73827E+03 

0.27068E+03 

4 

0.20992E-02 

0.19914E+04 

0. 59681 E+03 

0.27068E+03 

5 

0.21036E-02 

0.19955E+04 

0.52150E+03 

0.27068E+03 

6 

0.21063E-02 

0. 19981 E+04 

0.47460E+03 

0.27068E+03 

Element D(Rho)/DT 

D(P)/DT 

D(Vm)/DT 

D(Ve)/DT 

1 

0.83474E-01 

0.80936E+05 

0.26340E+06 

0.18311E+01 

2 

0.14618E-01 

0.14168E+05 

0.28578E+05 

0.18311E+01 

3 

0.57034E-02 

0.54621E+04 

0.80213E+04 

0.18311E+01 

4 

0.23664E-02 

0 . 22534E+04 

0.26043 E+04 

0.18311E+01 

5 

0. 10866E-02 

0.10325E+04 

0.98203 E+03 

0.18311E+01 


NOTES s 

Inner jet expansion, Thetal 0.09286 radians 
Outer jet expansion, Theta2 0.13457 radians 
Initial Jet with, BO, 0.01681 Ft 
Allowable Jet width, Bmax, 0.39000 Ft 

PREDICTED THRUST AND PERFORMANCE: 


Total ejector thrust = 
Isentroplc primary thrust = 
Thrust augmentation ratio = 


919.07166 Lb 
659.06592 Lb 
1.39451 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 

***** SINGLE POINT CALCULATION ***** 

T * 0.00070 sec 


FIELD VARIABLE PROFILE: 

Station -> Infinity IP IS 3 4 

AREA 11.70000 0.41998 9.32378 9.74376 17.55000 FT 

V 0.00000 860.78772 280.10501 301.62244 163.46561 FT/S 

P 2067.84009 2250.91992 1981.97961 1998.23413 2066.87256 LBF/FT2 

PO 2067.84009 3027.15869 2067.84009 2095.70337 2095.70337 LBF/FT2 

T 542.79999 698.86292 536.26273 553.28790 558.64166 DEG-R 

TO 542.79999 760.59998 542.79999 560.86810 560.86810 DEG-R 

RHO 0.07149 0.06044 0.06936 0.06777 0.06943 LBM/FT3 

MOOT 0.00000 21.85001 181.13047 199.17854 199.17854 LBM/S 

NOTES: 

Nozzle Pressure Ratio Is 1.46392 (Choked NPR Is 1.89293 ) 

FINITE VOLUME PROFILE 


Station Z 

B, Inner 

B, outer 

B.Jsl 

2 

0.18000E+00 

0 . 00000E+00 

0.41214E-01 

0.11696E-01 

3 

0.36000E+00 

O.OOOOOE+OO 

0.65618E-01 

0.15927E-01 

4 

0. 54000 E+00 

O.OOOOOE+OO 

0.90023E-01 

0.18886E-01 

5 

0.72000E+00 

O.OOOOOE+OO 

0.11443E+00 

0.21151E— 01 

6 

0.90000 E+00 

O.OOOOOE+OO 

0.13883E+00 

0.22978E-01 

Station Rho 

P 

Vm 

Ve 

1 

0.18785E-02 

0. 22509 E+04 

0.86079E+03 

0.28011E+03 

2 

0.20775E-02 

0.19712E+04 

0.11092E+04 

0.26850E+03 

3 

0.20937E-02 

0.19863E+04 

0.74354E+03 

0.26858E+03 

4 

0.21005E-02 

0.19926E+04 

0.59852E+03 

0.26859E+03 

5 

0.21041E-02 

0.19960E+04 

0.52205E+03 

0.26860E+03 

6 

0.21064E-02 

0.19982E+04 

0.47475E+03 

0.26860E+03 

Element D(Rho)/DT 

D(P)/DT 

D(Vm)/DT 

D(Ve)/DT 

1 

0.56055E-01 

0.55826E+05 

0.17555E+06 

-0.20917E+05 

2 

0.27180E-01 

0. 26487 E+05 

0.52730E+05 

-0.20917E+05 

3 

0.12250E-01 

0.11782E+05 

0.17126E+05 

-0.20917E+05 

4 

0.50163E-02 

0.47934E+04 

0.54920 E+04 

-0.20917E+05 

5 

0.17346E-02 

0.16524E+04 

0.15618E+04 

-0.20917E+05 


NOTES: 

Inner jet expansion, Thetal 0.09330 radians 
Outer jet expansion, Theta2 0.13476 radians 
Initial Jet with, BO, 0.01681 Ft 
Allowable Jet width, Bmax, 0.39000 Ft 

PREDICTED THRUST AND PERFORMANCE: 

Total ejector thrust = 909.19916 Lb 

Isentroplc primary thrust = 659.06592 Lb 

Thrust augmentation ratio = 1.37953 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 


***** SINGLE POINT CALCULATION ***** 
T = 0.00080 sec 


FIELD VARIABLE PROFILE: 


Station 

-> Infinity 

IP 

IS 

3 

4 


AREA 

11.70000 

0.41998 

9.32378 

9.74376 

17.55000 

FT 

V 

0.00000 

860.78772 

282.21085 

303.42273 

164.39374 

FT/S 

P 

2067.84009 

2250.91992 

1980.70361 

1998.74744 

2068.25391 

LBF/FT2 

PO 

2067.84009 

3027.15869 

2067.84009 

2097.42896 

2097.42896 

LBF/FT2 

T 

542.79999 

698.86292 

536.16406 

553.28845 

558.70764 

DEG-R 

TO 

542.79999 

760.59998 

542.79999 

560.95941 

560.95941 

DEG-R 

RHO 

0.07149 

0.06044 

0.06932 

0.06779 

0.06947 

LBM/FT3 

MOOT 

0.00000 

21.85001 

182.40828 

200.41866 

200.41867 

LBM/S 


NOTES: 

Nozzle Pressure 
FINITE VOLUME PROFILE 
Station Z 

Ratio Is 1. 
B, Inner 

46392 (Choked 
B, outer 

NPR Is 1.89293 
B.Jsl 

2 

0.18000E+00 

0.47418E-04 

0.41180E-01 

0.11992E-01 

3 

0.36000E+00 

O.OOOOOE+OO 

0.65550E-01 

0 . 16289E-01 

4 

0.54000E+00 

O.OOOOOE+OO 

0.89920E-01 

0.19109E-01 

5 

0.72000E+00 

0.00000E+00 

0.11429E+00 

0.21208E-01 

6 

0.90000E+00 

O.OOOOOE+OO 

0 . 13866E+00 

0.22894E-01 

Station Rho 

P 

Vm 

Ve 

1 

0.18785E-02 

0.22509E+04 

0.86079E+03 

0.28221E+03 

2 

0.20811E-02 

0.19748E+04 

0.11206E+04 

0.27057E+03 

3 

0. 20969 E-02 

0.19894E+04 

0.74972E+03 

0.27066E+03 

4 

0.21024E-02 

0.19945E+04 

0.60126E+03 

0.27069E+03 

5 

0.21051E-02 

0.19970E+04 

0. 52321 E+03 

0.27069E+03 

6 

0.21070E-02 

0. 19987 E+04 

0.47524E+03 

0.27070E+03 

Element D(Rho)/DT 

D(P)/DT 

D(Vm)/DT 

D(Ve)/DT 

1 

0.36388E-01 

0.35501E+05 

0.11427E+06 

0.21058E+05 

2 

0.31656E-01 

0.30789E+05 

0.61846E+05 

0.21058E+05 

3 

0.19465E-01 

0.18730E+05 

0.27408E+05 

0.21058E+05 

4 

0.10562E-01 

0. 10093 E+05 

0.11643E+05 

0.21058E+05 

5 

0.53912E-02 

0.51329E+04 

0.48831E+04 

0.21058E+05 


NOTES: 

Inner Jet expansion, Thetal 0.09286 radians 
Outer jet expansion, Theta2 0.13457 radians 
Initial Jet with, BO, 0.01681 Ft 
Allowable Jet width, Bmax, 0.39000 Ft 

PREDICTED THRUST AND PERFORMANCE: 


Total ejector thrust * 
Isentroplc primary thrust * 
Thrust augmentation ratio = 


920.05444 Lb 
659.06592 Lb 
1.39600 
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APPENDIX D - OUTPUT LISTING FOR CASE 1 


***** SOLUTION PROFILE ***** 


T 

DMpDT 

THRUST 

(o-sec) 

(Ibm/s) 

(lb) 

0.40000 

18.70000 

926.57001 

0.50000 

18.70000 

926.50159 

0.60000 

21.85001 

919.07166 

0.70000 

21.85001 

909.19916 

0.80000 

21.85001 

920.05444 


PHI 


1.82127 

1.82113 

1.39451 

1.37953 

1.39600 
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TABLE I. - EJECTOR GEOMETRY FOR STOVL TEST CASE 
[Configuration for 12 nozzles; 1 12-in. wide duct; 3 orifices per nozzle.] 


Station 

Direction 

X 

Y 

Z 

Characteristic dimension, ft. 

0 

9.333 

1.875 

0.900 

1 

.108 

.108 

.900 

2 

.778 

.347 

- 

3 

9.333 

1.042 

2.917 

4 

9.333 

1.875 

- 


Mixing region 


Volume, dV 

Area, dA 

Scale factors 



Number of Units 

For dV and dA 



n 

n 

N = n n 



X 

y 

x y 

Wdxdz = Zj 

Wdx 

12 

3 

36 

Wdx = ZjA 
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Figure 1. - General program structure. 
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Figure 2 - Outline of the SPC solution procedure. 















VtewA-A 



Figure 3. - Geometric ejector approximation. 
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Normalized ejector response (magnitude Is dB) 




Feed 


Figure 4. - STOVL ejector application. 



Figure 5. - Quasi-steady response criteria. 
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System thrust (XI 000 lb) 



Elapsed time (milliseconds) 

Figure 6. - Result from transient flow test case. 
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